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Abstract 

Quarks and Antiquaries in Nuclei 
Jason Smith 

Chair of Supervisory Committee: 

Professor Gerald Miller 
Department of Physics 

The Chiral Quark-Soliton model provides the quark and antiquark substructure of the nu- 
cleon, which is embedded in nuclear matter. This provides a new way to asses the effects 
of the nuclear medium on the nucleon. We calculate nuclear binding and saturation, de- 
scribe the European Muon Collaboration effect consistently with Drell-Yan experiments, 
and predict modifications to the polarized parton distributions. We also calculate the elec- 
tromagnetic form factors of a bound proton, and find significant modifications of the elec- 
tromagnetic form factors in the nuclear environment. In every case, the properties of the 
sea of quark-antiquark pairs in the nucleon are very important, and serve to mitigate the 
valence quark effect. The changes in the sea quarks when the nucleon is immersed in the nu- 
clear medium are the primary mode by which consistency is maintained with experimental 
constraints (Drell-Yan data, magnetic moment), while still maintaining a significant effect 
needed to explain the deep inelastic scattering and polarization transfer data. 
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Chapter 1 
INTRODUCTION 



The past century has seen the birth and incredible progress of the physics of hadrons 
and their strong interactions from the discovery of the nucleus |88| to the recognition of 
Quantum Chromodynamics (QCD) as the fundamental theory of the strong interactions 
with the 2004 Nobel prize in physics |54l I84j . Despite being the fundamental theory, de- 
scribing a proton, much less an entire nucleus, in the context of the quark and gluon degrees 
of freedom of QCD is still out of reach. The fundamental theory is non-perturbative at the 
low energy, long wavelength scales of nuclear physics, but it does tell us some very useful 
information about the physics of hadrons. It gives us a scale Kqcd that separates pertur- 
bative and non-perturbative physics, tells us how to translate our knowledge at one scale 
to another through the evolution equations, gives us an approximate chiral symmetry im- 
plying its spontaneous breaking at low scales, and hints at the non-perturbative degrees 
of freedom with the instanton solutions to the QCD field equations providing a possible 
mechanism for that symmetry breaking |104j . The enormous effect of chiral symmetry has 
been encapsulated in the very successful low energy effective theory, Chiral Perturbation 
Theory (xPT) 

These high and low energy theories have yet to meet, and a frontier of hadronic physics 
lies in the intermediate range of length scales available to present experiments where neither 
QCD nor xPT have useful perturbative expansions. Neither fundamental quarks and glu- 
ons nor point-like hadrons provide a complete description, so including the non-perturbative 
information that hadrons are bound states of valence quarks in a polarized vacuum is neces- 
sary. One way to probe these intermediate length scales and this non-perturbative physics 
is to examine the short distance structure of a large object. The prime example is the 



2 



European Muon Collaboration (EMC) effect measured in lepton-nucleus Deep Inelastic 
Scattering (DIS) where the short distance (~ 5 GeV, or ~ 10~ 2 fm) structure of nuclei 
differs from that of a collection of free nucleons. This measurement showed that the quarks 
and antiquarks in bound nucleons carry less momentum than in free ones, and implied 
that the medium modification of nucleon structure could be significant for many nuclear 
observables |5()j . Explanations of the EMC effect have remained a diverse and controversial 
theoretical subject. Until a measurement of dimuon production in proton- nucleus scattering 
through the Drell-Yan process [5] demonstrated separately that the antiquarks in nucleons 
are virtually unchanged in a nuclear system, there was virtually no other handle on medium 
modifications in DIS. Many explanations either fail to describe the Drell-Yan data, fail to 
address antiquarks at the non-perturbative scale of nuclei, or fail to address antiquarks 
altogether. 

There is a hint to the solution from evolution in QCD; dynamical rescaling 29 suggests 
that the EMC effect is due to a lowering of the starting scale for the renormalization group 
equations. While successful at describing the phenomenology of both the EMC effect and 
Drell-Yan experiments, like the situation for QCD and chiral symmetry breaking, it is silent 
on the actual mechanism and scale. 

The first evidence that the protons and neutrons had structure came from the measure- 
ment of their electromagnetic form factors in elastic lepton- nucleon scattering [Jl] . Recent 
polarization transfer experiments at TJNAF (Thomas Jefferson National Accelerator Facil- 
ity, also referred to as JLab) jlU2j observed a difference in the electromagnetic form factors 
of a proton bound in a Helium nucleus compared to a free one. This seems to add to the 
evidence in favor of the modification of hadrons in the nuclear medium. 

We would like to address the subject of medium modifications with a theory that is 
consistent with QCD and xPT. Any description of the EMC effect must be consistent with 
the constraints set by both DIS and Drell-Yan data. Thus a successful model must include 
antiquarks as well as quarks, and show how the medium modifies both the valence and sea 
quark distributions. Our purpose is to provide a mechanism for that modification within the 
Chiral Quark-Soliton (xQS) model |681 1161 EM I26j . It is essentially a relativistic mean field 
approximation approach to baryons motivated by the Instanton Liquid Model (ILM) of the 
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QCD vacuum and the large Nc expansion (discussed in Chapter^. This phenomenological 
model has many desirable qualities: the ability to describe a wide class of hadron observ- 
ables with surprising accuracy given its simplicity, the inclusion of antiquarks, the positivity 
of quark distributions [HE]) polynomiality of Generalized Parton Distributions (GPD) |94j . 
and a basis in QCD 39 . The model also predicted the mass and narrow width jlj of 
the 9 + pentaquark for which there was encouraging early evidence [BIE^, but remains a 
controversial subject with many experiments showing null results |43j . One of the draw- 
backs of the model is that, in the form presented here, it lacks explicit confinement. This 
is the property of QCD that quarks and gluons are never seen as free particles. The lack of 
confinement in the xQS model has not prevented it from describing a wide variety of free 
nucleon observables Here we show how the model describes nuclear saturation prop- 
erties, reproduces the EMC effect, satisfies the bounds on nuclear antiquark enhancement 
provided by Drell-Yan experiments, and predicts large effects in spin polarized structure 
functions and electromagnetic form factors. 

We will continue this introduction to the subject with a more detailed description of deep 
inelastic scattering and the EMC effect in Section 11.11 and dimuon production in Section 
11.21 We will then present the xQS model in Chapter^ and use it to build a nuclei out of 
quarks and antiquarks in Chapter |3 The medium modifications of the quark distributions 
and electromagnetic form factors are calculated in Chapters 0] and [51 respectively. In the 
Appendices, we elaborate on the connection of the xQS model to xPT and QCD via lattice 
calculations of the non-perturbative quark propagator, provide background on the small 
binding effect predicted by relativistic mean field theory for finite nuclei and nuclear matter, 
and present some details for the numerical work. This work is built on results that have 
been published previously elsewhere [211 EHl EE3 [100 . 

1.1 Deep Inelastic Scattering and the EMC Effect 

Deep inelastic lepton-nucleon scattering occurs when a virtual photon, mediating the inter- 
action, of high energy and momentum q (—q 2 3> ^qcd) fragments the target nucleon N 
of mass Mjv into a host of particles labelled X. The electromagnetic interaction couples 
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Figure 1.1: Basic diagram of the deep inelastic scattering process fx + N —* fi + X with 
momentum transfer q. 

to all of the electrically charged constituents in the target nucleon. Deep inelastic scatter- 
ing experiments see the valence quarks as well as the virtual quarks and antiquarks in the 
polarized QCD vacuum simultaneously. 

The case where the lepton is a muon is illustrated in the Feynman diagram in Fig. 11.11 
The use of the optical theorem allows us to package our ignorance of the photon-nucleon 
vertex and fragmentation function in Fig. 11.11 In inclusive inelastic scattering (where one 
measures the cross section with any final state X), we can write the sum on the final states 
X of the diagram in Fig. 11.11 as 




= 2 Im W» v {k,q) 



We can use Lorentz covariance and current conservation to write down the form of W^ v 

W^(k,q)=^-^^W 1 (k,q) + ^^-q^^ (V - q v ^j W 2 {k, q), (1.1) 

but we need a model to calculate , specifically the functions W\ and W%, explicitly. 
For example, we can use the parton model of Bjorken and Feynman as in Fig. 11.21 where 
the nucleon is composed of massless (or nearly massless) point-like, non-interacting, mostly 
spin-1/2 constituents. This model is justified by the asymptotic freedom |84[ I54j of QCD, 
where at high energy the strong interaction coupling a s between the quarks and gluons goes 
to zero. 
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Figure 1.2: The handbag diagram in the parton model. The virtual photon 7* transfers 
momentum q to the nucleon. The nucleon has momentum k, and the struck quark has 
momentum p. 



It is conventional to use the momenta in Fig. 11.21 to define the variables 

Q 2 = ~q 2 (1.2) 

k ■ q , 

- -« 2 « 2 -* + (1.4) 



2k -q 2M N v k+ 

where x is called the Bjorken variable. It has the interpretation in the parton model of 
the fraction of the nucleon momentum carried by the struck parton. The parton model is 
defined by the constituents with momentum p in Fig. 11.21 satisfying 

(p + qf = p 2 (1.5) 
=>2p-q-Q 2 = (1.6) 

so that if p = xk, then x is given by Eq. (|1.4I) . We have also given x in terms of light 
cone variables, A = A ± A± = (A 1 , A 2 ) that are closely connected to the infinite 
momentum frame where the parton model is defined. The approximation in Eq. 1)1. 4J) is 
that at large Q 2 and v, the struck quark has effectively zero transverse momentum after 
the interaction with the photon. 

The Bjorken limit ^7] is defined as Q 2 , v — > 00 with Q 2 fv finite, whereupon the functions 
Wi exhibit Bjorken-Feynman scaling, and we define the structure function of the nucleon 



Mn w2(v,q 2 ) ^Y.^m^ qi{Q2,2MNV) - F2{N){x) (L7) 
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where i runs over the quark flavors. This scaling breaks down experimentally as the function 
acquires additional Q 2 dependence. The Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) evolution equations [121 ESI El (see Appendix EJ) derived from QCD describe 
the scaling violation by allowing one to calculate how the distribution of momentum qi(x) 
among the quarks (and gluons) changes with the momentum scale Q 2 . It is beyond the 
parton model, and is one of the most successful results of QCD. The structure function 
Eq. I)1.7j) is defined in terms of the quark distribution functions qi(x), which, in the parton 
model in Fig. 11.21 is given by 

qi(x) = k+ J ^ e^^'ikm^^z-^k) (1.8) 

The tp are quark field operators, and ip = 7/^7° (we will use ^ for nucleon field operators). 
There are two interpretations of quark distributions: 

• Probability distributions of partons carrying fraction x of the nucleon momentum in 
the infinite momentum frame, and 

• QCD matrix elements of quark bilinears with light-like separation and the Operator 
Product Expansion (OPE) for the moments of the distributions. 

We have so far used the first one, but these two definitions are equivalent in the xQS model 
because Feynman's hypothesis that the transverse momenta of partons do not grow with 
Q 2 is satisfied 36 . 

The second approach is advantageous for the discussion of dynamical rescaling |29| . We 
can define the moments of the structure function of a nucleus F^ 1 with momentum P in 
Wilson's operator product expansion 

Mi A \Q 2 ) = ^j\xx- 2 F^(x,Q 2 ) 

= ^£^WVi)<5$(^) (1-9) 

n,i 

where the sum runs on all twist- 2 operators contributing to F% (twist of an operator is given 
by t = di — n, where dj is the mass dimension of the operator), and O is the matrix element 
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(p\o\Z..»Ma)\ p ) = ofliridP* ■ ■ ■ p ^n (i-io) 

The essence of the dynamical rescaling description of the EMC effect is that the renormal- 
ization scale for nuclei is less than that for a free nucleon 

fJ-A < Mi < Q 2 

so that the nuclear operators 0^ A ' undergo 'extra evolution' going up to the DIS scale Q 2 , 
which accounts for the EMC effect. QCD does not tell us what ii\ is, nor its relationship 
to but numerical studies {22J show that this process can describe both the EMC effect 
and Drell-Yan experiments |2j. 

The Wilson coefficients c in Eq. in the OPE are dependent on the starting renor- 
malization scale, and the essence of the dynamical rescaling picture of the EMC effect is 
that the scale is different in nuclei than for a free nucleon. Nuclear interactions are effective 
QCD interactions which renormalize a free nucleon into a bound one, but as we cannot 
yet in QCD calculate the parton distributions of a free nucleon, much less one in a nuclear 
medium. We need a model to probe the mechanism of that 'renormalization'. 

In the handbag diagram Fig. 11.21 the struck quark in the intermediate state is renormal- 
ized by QCD from the same starting scale as the other quarks in the nucleon, and hence 
all feel the effects of nuclear interactions at the low scale. In Chapter 0] we will calculate 
the effect of the presence of other nucleons on the quark distributions. One can imagine 
that we are calculating the effect of different starting scales for free and bound nucleons by 
showing what a bound nucleon parton distribution looks like at the scale of a free nucleon. 

The use of a manifestly covariant formulation to derive the expression for the nuclear 
structure function appearing in Eq. (jl.9jl leads to a convolution formula jSH] 

r 2 



- = / dyf(y)Fi N> (x/y) 



A , _ 2 ,-,„ (LID 

The function f(y) is the nucleon momentum distribution in the nucleus, and is given by 
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where the function X P) 1S the Fourier transform of the connected part of the nucleon 
Green's function. The calculation of Eq. (j!.12|) is detailed in Appendix IU1 

The properly normalized nuclear structure function F2 (x)/A is smaller than the free 
nucleon structure function F% (x) for values of 0.3 < x < 0.7. This is the EMC effect 0, 
and it has been known for more than twenty years. The underlying mechanism responsible 
for the change of the distribution of momenta among the constituents of the nucleus remains 
beyond our ability to calculate in QCD. We must resort to models, and there is a clear 
interpretation within the parton model: a valence quark in a bound nucleon carries less 
momentum than a valence quark in a free one. One popular mechanism for this reduction 
in momentum involves ordinary nuclear binding which, in its simplest form, is represented by 
using the nucleon momentum distribution Eq. (|1.12|) to evaluate the free nucleon structure 
function at a value of x increased by a factor of the average separation energy divided by 
the nucleon mass e/M/y — 0.04. This picture encapsulates any nuclear effects in the nuclear 
structure F% /A in the function f(y). 

The validity of the binding effect has been questioned in the reviews [HI 023 ESI HE]; 
and we have shown |79[ I98j that conventional relativistic nuclear mean field theory fails to 
describe the phenomenology of the EMC effect in terms of f(y) in Eq. (jl.llj) . These results 
are summarized in Appendix [DJ We believe that this failure to describe the EMC effect 
provides additional motivation to study quarks and antiquarks in nuclei beyond trying to 
connect QCD to nuclear physics, and to look for medium modifications of the bound nucleon 
in the only other part of Eq. Ql.ll|l available: the structure function F^ N ^ ■ We must also 
remain consistent with the Drell-Yan experiments, which we describe next. 

1.2 The Drell-Yan Process and Dimuon Production 

The proton-nucleus Drell-Yan cross-section can tell us about the antiquark content of the 
nucleus, and hence serves as another handle on nuclear parton distributions. The exper- 
imental data are consistent with no modification of the antiquark distribution [2j. Many 
models fail to describe, or even address, dimuon production in the Drell-Yan process de- 
picted in Fig. 11.31 Dynamical rescaling matches the trend of the data [2], but, as in its 
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description of DIS data, it does not give any insight as to the mechanism. If x p is the 
momentum fraction carried by the (anti)quark in the proton, and xa is the same for the 
nucleus, then one finds the cross section from Fig. 11.31 12 

d 2 a A7ra 2 K 1 



E e * M P) (* P , Q 2 )#W Q 2 ) + tf\x P , Q 2 )i A \x A , Q 2 ) 



dx p dxA 9s XpXA 

(1.13) 

where a ~ 1/137 is the electromagnetic coupling, and y/s is the energy in the center-of-mass 
frame of the proton and nucleus. The factor K is an enhancement factor accounting for 
other dilepton production modes independent of x p and xa and the target and projectile, 
see |12j . and references therein. If we look experimentally at x p > xa + 0.2, then at high 
Q 2 the first term dominates (the antiquark distribution for the proton is effectively zero 
for x p > 0.2 and Q 2 large). This can be achieved by controlling the kinematics of the 
experiment, so we can use the Drell-Yan process to examine the antiquark content of the 
nucleus. In order describe the experimental data, we must have a model that includes both 
quarks and antiquarks at the model scale 1 . We introduce such a model in Chapter [2 

The evolution equations in Appendix [H] generate antiquark, as well as gluon, distributions as the scale 
changes. While different starting distributions will generate different antiquark distributions, a lack of 
antiquarks in the model ignores possible nuclear effects at the starting scale. A common reason for failing 
to describe the Drell-Yan data is that models produce too large an enhancement in the antiquark content, 
so addressing antiquarks with a model at the model scale is essential. 
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Chapter 2 

THE CHIRAL QUARK-SOLITON MODEL 

We will begin our study in earnest by providing motivation for the xQS model, following 
the reviews |39[ I34| . First, we will discuss instantons, and how they break chiral symmetry. 
Constructing the model itself, and deriving some of its basic properties in the meson and 
baryon sectors is next. The meson sector connects us to xPT, and the baryon sector connects 
us to our description of nuclear elastic, and deep inelastic scattering experiments. 

2.1 Results From The Instanton Vacuum 

Quantum Chromodynamics is a Yang-Mills gauge theory with gauge symmetry SU(Nc), 
where the number of colors is Nc = 3. In Euclidean space 1 , the QCD field equation for the 
gluons contains solutions called instantons. These solutions never appear at any order in 
the perturbation theory expansion in the strong coupling ot s , and hence may provide insight 
into non-perturbative QCD. Instantons can exist concentrated at a single point in spacetime 
and are sometimes called pseudoparticles as they do not follow spacetime worldlines like 
ordinary particles. Lattice simulations provide evidence that the QCD vacuum is populated 
by instantons with a size p ~ 0.3 fm, which is approximately the scale of chiral symmetry 
breaking, and a separation R = (N/V) 1 ^ ~ 1 f m where N/V is the instanton number 
density. The instanton packing fraction, equal to the four dimensional instanton volume 
7r 2 /9 4 times the instanton number density, is then 

ttV _ 1 
i? 4 ~ 8' 

which gives a small parameter aside from a s in perturbative QCD to use in expansions. This 
is the basis for the Instanton Liquid Model (ILM) picture of the QCD vacuum. Though 

1 Euclidean spacetime related to ordinary Minkowski spacetime with time t by taking t — ir with r being 
the Euclidean time. 
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dilute based on the packing fraction, the instantons are still strongly interacting, therefore 
a liquid, not a gas, picture is more accurate. 

Because of their topological properties, the instantons are accompanied by 2Nf zero 
modes of chiral fermions (Nf each of left and right handed fermions), where Nf is the 
number of flavors of quarks. Instantons are localized near a single point in spacetime and 
have 2Nf fermions meeting at that point, so they look like interaction vertices appearing in 
Feynman rules. In the case of Nf = 2, we have effective four-quark interactions. These are 
the so-called 't Hooft vertices |ll)4| ; a Feynman diagram for one with Nj = 2 (which we will 
use throughout) is in Fig. 12.11 A random ensemble of instantons jSZ] provides a mechanism 
to break chiral symmetry via these zero modes. Effectively one creates a helicity- flipping 
term in the Lagrangian that behaves like a mass term via scattering of chiral quarks by the 
't Hooft vertices. They create a dynamical quark mass 



M(p) = \J ?^<j>(pp/ 2 f (2.1) 

<P(z) = 2z I (z)K l {z)-I l {z)K (z)-- z h(z)K l (z) (2.2) 

where <j)(z) is the Fourier transform of the quark zero mode wave function bound to each 
instanton, N/V is the instanton density, and Nq is the number of colors in the gauge theory. 
The constant A is obtained from a detailed numerical calculation |34| . which results in a 
dynamical constituent quark mass with M(0) ~ 350 MeV. The dynamical quark mass falls 
to zero for p> 1/p ~ 650 MeV which acts as an effective momentum cutoff for the theory. 

If one combines a quark and an antiquark line from a 't Hooft vertex Fig. 12.11 one obtains 
a quark-meson vertex pictured in Fig. l2.2l where the type of meson depends on the quantum 
numbers. If one considers the correlation function of two quarks 

I?(p) = -J d 4 xe~ ip - x tr SOr,0)rS(0,x)r (2.3) 

T = 1,75,7 M ,757 M ,<V ( 2 - 4 ) 

where S(x, y) is the quark Green's function and the bar indicates an average over the 
instanton ensemble in the ILM, one finds a pole in the pseudoscalar channel at p 2 = cor- 
responding to our massless pions resulting from the dynamical breaking of chiral symmetry, 
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Figure 2.1: A 't Hooft vertex |104| for Nf = 2. The instanton is labelled by /, and the 
quark zero mode wave functions are represented by form factors \f~M. The arrows indicate 
the helicity of the quarks. 



and in the scalar channel at p 2 ~ 4M 2 (0) which is taken as the mass a scalar meson. There 
are no poles in the vector or tensor channels pUj. We will denote scalar, pseudoscalar, and 
vector with subscripts s, ps, and v in the following. 

2.2 The Model Lagrangian 

Using the quark meson vertex in Fig. 12.21 we can write down a Lagrangian with con- 
stituent (anti)quarks ip, ifi and pions tv(x) 

C$,M = ^ipiP - j d i p 1 d i p 2 e i ^-P 2 > x ^( Pl ) v ^M{^U^ 

W*(x) = exp^ 75 ^y^) (2.5) 

where f n ~ 93 MeV is the pion decay constant and Tj = <7j/2 are the isospin operators (<Tj are 
the Pauli spin matrices). This theory exhibits asymptotic freedom like QCD; the interaction 
disappears at large momenta with the function M(p) falling to zero leaving a theory of free 
massless quarks. We would like to replace this complicated nonlocal theory with a local 
one; this is done by replacing the momentum dependent mass M{p) with a constant M. 
The ultraviolet cutoff M (p — > oo) — > is replaced by a Pauli- Villars regularization fixed by 
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Figure 2.2: Quark-meson vertex formed from a 't Hooft vertex (Fig. 12.11) . where T = 
1) 75) 7^! 757^1 represent the channels onto which the quark zero modes are projected. 



the physical value of 




fl = 4iV c j 



d 4 k M 2 M 2 M 2 V 



(2vr) 4 (k 2 - M 2 ) 2 M 2 (k 2 - M 2 ) 2 



NpM 2 m 2 v 
= ~4ir 2 ~ ~M 2 ~ (2 ' 6) 

Using M = 420 MeV, the resulting value for Mpy ~ 580 MeV is comparable to the instan- 
ton size cutoff, and the resulting mass function is close not only to the instanton model 
Eq. (|2.2|) . but to an extrapolation of quenched lattice data to the chiral limit (see Appendix 
|A")) . The three functions are plotted in Fig. 12.31 This value for the constituent quark mass 
has also been shown to reproduce many of the properties of free nucleons |26j . which will 
be discussed in Chapters |1] and |SJ We introduce the notation 

M 2 

A-PV = A{M) - —^A(Mpv) 

as a shorthand to represent the Pauli-Villars subtraction, where A is any logarithmically 
divergent quantity in the xQS model. We use this method to regulate the divergences as it 
preserves the completeness of the quark states, which is important in the calculation of the 
quark distributions [37)] . 
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Figure 2.3: The mass function M(p) using one Pauli-Villars regulator fit to the pion decay 
constant (heavy line), the instanton model (light line) jSZ] and the extrapolation of fits to 
quenched lattice data |86| to the chiral limit using three Pauli-Villars propagators (heavy 
dashed line), see Appendix El 
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In the case of the nucleon state the pion field is taken to be the 'hedgehog' ansatze with 
profile function 0(r) and unit winding number. 

U(x) = exp(m-re(r)) (2.7) 
6(0) = -7T (2.8) 
0( r _> oo) = o (2.9) 

Unit winding number is guaranteed by the condition Eq. I|2.8j) . The local version of the 
theory ()2.5|) with constant mass, Pauli-Villars regulator and hedgehog profile function results 
in the Lagrangian for the Chiral Quark-Soliton (xQS) model 

C = ^{ip -Me !75 "- Te(r V (2.10) 

The pion field configuration corresponds corresponds to the soliton (Skyrmion) in the 
Skyrme model, hence its appearance in the name of the model (quarks, however, do not 
appear as degrees of freedom in the Skyrme model). This is a non-renormalizable effective 
field theory, so that the ultraviolet cutoff, the Pauli-Villars mass, has a physical meaning 
(observables can depend on it), namely the the instanton size. The constituent quark mass 
M could be taken to be a free parameter when we tackle medium modifications in Chapters^ 
and|SJ but we have chosen to fix it at 420 MeV. This value best reproduces a wide variety of 
nucleon observables [271 126j . and produces a mass function similar to Eq. (|2,2j) in Fig. 12.31 
It is worth noting that we could have just written down the Lagrangian (|2.10j) without 
motivation from the instanton picture of the QCD vacuum as the simplest field theory of 
constituent quarks consistent with chiral symmetry (a chiral transformation if) — > e ll5a ^ip 
can be absorbed into a redefinition of the pion field 7r). 

Quarks are described by single particle wave functions ip n that are solutions to the Dirac 
equation 

{ip - Me' 75,lTe(r) )i(r)e- iE "' = (2.11) 
e ,E " t t(r)(i^-Me i75,lTe(r) ) = (2.12) 

with energy eigenvalues E n . They are represented as a discrete set labelled by n for sim- 
plicity, but the spectrum of the Dirac operator contains a single bound state (which we will 
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refer to as the valence level, and label with a superscript v), as well as a positive and nega- 
tive energy continua for energies \E\ > M. At most one bound state is guaranteed the the 
unit winding number Eq. (|'2.8j) of the hedgehog field (there may be no bound state solution 
for certain values of M). Unit winding number is a necessary, not sufficient, condition to 
obtain unit baryon number given by the number of quarks in the valence level; the winding 
number is not the baryon number as it is in the Skyrme model. 

The states in the negative continuum are filled in both the vacuum sector with U = 1 
and the nucleon sector with profile function 0(r). For any observable, one must subtract 
the contribution from the vacuum (U = 1) to obtain a physical result. These observables 
can still be logarithmically divergent after the vacuum subtraction, requiring Pauli-Villars 
regularization. It is the filled negative energy continuum that allows one to describe the 
nucleon quark sea: the virtual quarks and antiquarks that appear out of the polarized QCD 
vacuum in the presence of a nucleon. Positive energy antiquarks appear as a hole: a lower 
density of states in the negative energy continuum. Conversely, negative energy quarks 
appear as a higher density of states. The pions in the theory are lighter than 2M because 
they are density fluctuations in the negative continuum consisting of a positive energy 
antiquark and a negative energy quark so that their mass is given by m n ~ M — M = 



There is an equivalence in calculating observables where one can use sums over the 
occupied states (the valence level and negative energy continuum) or the unoccupied states 
(the positive continuum) [SOI- This equivalence is deeply rooted in the fact that the model is 
a relativistic quantum field theory, and is needed to preserve general conditions such as the 
positivity of the antiquark distribution and the baryon and momentum sum rules discussed 
in Chapter 0J It can be spoiled by the regulator, such as Schwinger's proper-time regulator 
or sharp cutoffs, but Pauli-Villars regularization preserves the completeness of the quark 
states and maintains the equivalence |36j . 

The field equation for the profile function follows directly from the Lagrangian (|2.1U|) 
and is given by 




(2.13) 
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where the quark scalar and pseudo-scalar densities are 



(2.14) 



n 



r = 



l,ry 5 n • r 



For a single nucleon this sum is over all occupied states which includes the negative energy 
continuum and the single valence level. In the nuclear medium, this becomes more com- 
plicated, and the general case would be exceedingly difficult. In symmetric infinite nuclear 
matter, with equal constant densities of protons and neutrons filling all space, we will be 
able to make some simplifications. This is an often used approximation to describe heavy 
nuclei QUE]. 

2.3 The Nucleon and the Large Nq Limit 

Aside from the gauge coupling used in perturbative QCD, there is another 'small' dimen- 
sionless parameter that we can expand in: the inverse number of colors Nq . It is believed 
that many of the important properties of QCD follow into the limit where Nq is large; we 
will leave the basics of this expansion regarding baryons to Ref. jlU9j . We will use the limit 
here to justify the xQS model as a relativistic mean field theory of baryons [SH]) simplifying 
the calculation of the nucleon mass, form factors and quark distributions. We will follow 
Ref. jHH] to obtain an expression for the nucleon mass. 

In Euclidean space (in which we will work in this section), the partition function for the 
xQS model is 



The integration over the quark fields can be done, giving the effective action for the pion 
fields (including the vacuum subtraction) 




(2.15) 




(2.16) 



One can expand this in derivatives 




(2.17) 
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and see that one obtains the 0(p 2 ) term in the momentum expansion in %PT (note this also 
defines Eq. (|2.6j) for f n ). If one includes a current quark mass, one obtains a Lagrangian 
consistent with the 0(p 4 ) Gasser-Leutwyler Lagrangian and large Nc- This is elaborated on 
in Appendix but we have already made the connection between xPT an d QCD through 
the xQS model. 

From the effective action, we can also determine the energy in the background pion field 
in the limit of long Euclidean times T. We write Eq. (|2.1(i|) in terms the Hamiltonian H, 
and with U = 1 



ScffM 



-Nc Tr 



log (id 4 + iH) - log [id A + iH^ 

log (E + iH) - log (e + iH^ 
E E 



' ' ' 2vr \E + iH E + iHW 



(2.18) 
(2.19) 

(2.20) 



KcTl ^2 E n - y: em 

\E n <0 
-NcTEaaAhr] 



4 0) <o 



(2.21) 
(2.22) 



One expects the limit of the correlation function of two nucleon currents to give the 
nucleon mass at long times T 

lim ( J(0, T) jt(0, 0)) ~ e~ MNT (2.23) 

T — >oo 

where the nucleon current is given by 

J(x) = ^ ai "^Si lA (z) • • • i> aNcfNc (x) (2.24) 

where on are color indices, fi are flavor indices, J and J3 are angular momentum quantum 
numbers, and T and T3 are the isospin. We subsume the description of the nucleon quantum 
numbers in terms of its constituent quarks into the matrix T. The correlation function is 
defined in terms of the partition function, and the integration over the quark fields can be 
done. This results in the effective action and the quark Green's function in the background 
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pion field 



(J(0,T)jt(0,0)) = \ \ DirDi>D^J(0,T)J^(0,0)e~ s ^^ 

N C 

^{f} r M / 



r {JT} r {jT } / ^nG /ift (T,0|vr)e- 



The long time limit of the quark Green's function in the background pion field gives us the 
valence energy level 



lim G /ift (T,0|vr)~e- T ^W (2 . 25) 

1 — >oo 



Putting this together we obtain 



lim <J(0,T)jf(0,0)) ~ Dne-TNciwW+EhMW) 

T — >oo 



which gives us, in the saddle point approximation (leading order in Nc), the nucleon mass 

M n [tt{x)] = min N c (£>] + E &eld [ir}) (2.26) 

{ttO)} 

We have neglected quark loop corrections at large Nc in the saddle point approximation (and 
therefore quantum fluctuations of the pion field, which propagate through quark loops). In 
this sense, the large Nc limit effectively represents a relativistic mean field approximation. 
The abstract functional Eq. (|2.26|) is minimized by solving the system self-consistently with 
the pion field replaced by the profile function 0(r). This is the subject of Chapter 01 
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Chapter 3 



QUARKS AND ANTIQUARKS IN NUCLEI 



In order to describe nuclei in terms of the quarks and antiquarks in the xQS model, 
there needs to be a way to connect the nuclear physics and the quark physics. This has 
been done in the context of the MIT bag model [HJ to create a so-called Quark-Meson 
Coupling (QMC) model. We will use some of the ideas of that model, but also include some 
phenomenology of medium modifications and physics of the QCD vacuum. We begin with 
the role of the vacuum quark scalar condensate, and then construct a two-tier self-consistent 
calculation. We are able to obtain nuclear binding and saturation, and reproduce standard 
properties of nuclear matter. We therefore create a model of nuclei that includes the physics 
of valence and sea quarks. 

3.1 The Effect of the Medium 

We will begin with some motivation for our procedure to couple the quark substructure of 
the nucleon to the nuclear medium. Through the use of QCD sum rules, Ioffe derived a 
relationship between the vacuum scalar condensate, (ipip)o, and the nucleon mass. One can 
re-derive this estimate in a constituent quark field theory such as we are using here. We 
begin with the scalar condensate 



where the divergent loop integral is regulated by a momentum cutoff (playing the role of 
the Borel mass in the QCD sum rule approach). Using the fact that constituent quarks are 
essentially defined as having a mass ~ Mn/Nc, we can rewrite Eq. (|3.1j) as 




1 



N C MA 2 
4tt 2 



(3.1) 



M N 



(3.2) 
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Although Eq. (|3.2|) is not a very accurate estimate, it does highlight the role of the conden- 
sate. It will be modified in the presence of other nucleons. 

The condensate at finite density, proportional to kp where kp is the Fermi momentum, 
can be written in terms of the nuclear scalar density pg and the nucleon sigma term ctjv 
(SU as 



W)fc F = (tp^)o ~ (V^O-^TToPf • (3.3) 

mitt 



We can then substitute Eq. ()3.3|) into Eq. (|3.2|) to obtain a schematic picture of the effect 
of the nuclear medium on the nucleon mass 



2 



4-7T — 

M N {k F ) ~ -_(^) 



i _ ° N n 



(3.4) 



We will effectively replace the combination of the vacuum condensate, pion mass, decay 
constant and the the sigma term in Eq. (|3.4|) with a dimensionless scalar coupling constant 
g s , which we take as a free parameter. The value of g s is determined from a fit to nuclear 
density and binding energy. 

Using this dependence of the nucleon mass on the nuclear medium as a guide, we in- 
corporate the medium dependence in the model by simply letting the quark scalar density 
in the field equation ()2.13j) contain a (constant) contribution arising from other nucleons 
present in symmetric nuclear matter. This models a scalar interaction via the exchange of 
multiple pairs of pions between nucleons. We take the scalar density to consist of three 
terms: 1) the constant condensate value {if^ip)o (in the vacuum or at large distances from 
a free nucleon), 2) the valence contribution p v s and 3) the contribution from the medium 
which takes the form of the convolution of the nucleon and valence quark scalar densities 
as in the QMC model [ZH] 

p« s (r) ~ (W) + pUr)+g s Jdr'p?(r-r')p v s (r') (3.5a) 
= (H>)o + P v .(r)+g,P?(k F ) (3.5b) 

We define an effective condensate in the presence of nuclear matter by 

{^) + g s P s N (k F ). (3.6) 
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We take the pseudoscalar density to have only the valence term p\ s ~ Pp S ; the two other 
contributions analogous to the first and third terms of Eq. H3.5|) vanish due to symmetries 
of the QCD vacuum and nuclear matter. These approximations to the densities neglect 
the precise form of the negative continuum wave functions. We will see that the result- 
ing free nucleon profile function has little discernible difference from a fully self-consistent 
treatment (see Fig. 13.51 in Section l3.3j) . demonstrating the excellence of this approxima- 
tion. We provide a schematic picture of the scalar density in a free and bound nucleon in 
Fig. 13.11 Incorporating this dependence on nuclear matter in the field equation Eq. ()2.13j) . 
the in-medium field equation for the profile function becomes 



(i.e. medium effects serve to bring us closer to chiral symmetry restoration). 

The (constant) nucleon scalar density in Eq. (|3.5j) is determined by solving the nuclear 
self-consistency equation |108| 



The dependence of the nucleon mass, and any other properties calculable in the model, 
on the Fermi momentum hp enters through Eq. (|3.8|) . Thus there are two coupled self- 
consistency equations: one for the profile, Eq. (|3.7|) . and one for the density, Eq. (|3.8|) . 
These are iterated until the change in the nucleon mass Eq. (|3.15|) is as small as desired (in 
our case, AMjy < 0.1 MeV) for each value of the Fermi momentum. 

We introduce a phenomenological vector meson field (with mass fixed at m v = 
0.77 GeV and coupling g v ) |108j exchanged between nucleons, but not quarks in the same 
nucleon (i.e. we ignore the spatial dependence of the vector field in the vicinity of a nucleon, 
treating only the nuclear mean field). In the mean field approximation (MFA) V l=1,2,5 = 0, 
and the vector potential has a simple field equation 108 





(3.8) 



V° = ^P v N (k F ) 



(3.9) 



where the vector meson couples to the vector density 




(3.10) 



23 




Figure 3.1: Cartoon of the scalar and pseudoscalar densities as a function of radius from 
the center of the nucleon. The inverse constituent quark mass sets the typical size scale. 
Pictured are (a) the vacuum condensate, (b) the effective condensate Eq. ()3.6|) . (c) the 
free nucleon scalar density, (d) the bound nucleon scalar density, and (e) the pseudoscalar 
density. 
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It is a proxy for uncalculated soliton-soliton interactions used to obtain the necessary short 
distance repulsion which stabilizes the nucleus. The resulting energy per nucleon is 

K 4 f' kp rfik , 1 n 2 

-A=-W-\\ 7^#+W+ ^(M. (3.H) 
A PB(kF) J (27iT 2 mf, 

The single-particle Dirac Hamiltonian for the quarks in the hedgehog pion field, after in- 
cluding the vector interaction, is 

H[@]4> n {r) = [-ij°~f k (d k + ig v V k )+g v Vo + M~f°(cosO + i~f 5 T ■ nsm@)]ip n (r) 

= E n ip n (r) (3.12) 

where the effect of the nuclear medium, the dependence on kp, enters through the vector 
potential, as well as the profile function © = 0(r, fcp). 

3.2 The Numerical Procedure 

The Dirac Hamiltonian is diagonalized in the discrete Kahana-Ripka basis |HZj (see Ap- 
pendix The basis functions fa are defined by the solutions to the free single-particle 
Dirac Hamiltonian 

H^fa(r) = H 7 V4 + M 7 ]fa(r) = Effair) (3.13) 
The full Hamiltonian is then solved by diagonalizing the matrix equation 

Y,H i3 [Q]c n] =E n c ni (3.14) 

3 

where 

H tJ [Q] = J dr4(r)H[e]^(r) 

and 

ipn(r) = ^2c ni fa(r) 

i 

The use of a finite basis turns the positive and negative energy continua into a discrete 
set of states. Our interpretation of these continua as quarks and antiquarks in the nucleon 
sea remains the same. Antiquarks (holes) appear as lower densities of discrete energy levels 
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in the negative energy states. In the finite basis, we can write the nucleon mass Eq. 1)2.26)1 

as 



Starting at zero Fermi momentum and an initial guess for the profile function as the first 
step in the self consistent calculation, the valence level wave function is used to calculate 
the quark scalar and pseudoscalar densities (there is no medium contribution at kp = 0). 
These are used in Eq. 1)3.7)) to calculate a new profile function. 

This new profile function then gives us new quark densities after diagonalizing the Hamil- 
tonian Eq. (|3.12j) . This procedure is iterated until self consistency is reached when the 
change in the valence energy level AE V < 0.01 MeV (the first self-consistent loop). Next, 
the eigenvalues are used to calculate Mjy. At kp = 0, the procedure is done here and kp is 
incremented a step Akp = 0.01 fm. 

With kp ^ 0, there is now a contribution to the scalar density from the medium in the 
profile field equation 1)3. 7|l . as well as a non-zero vector potential. We calculate the nucleon 
scalar density via Eq. 1)3.8)) . using the nucleon mass from the previous step, M^(kp — Akp), 
as the initial guess. The Hamiltonian is again diagonalized and the profile is calculated 
self consistently using the valence level and the medium contribution. The nucleon mass is 
then calculated, and the nucleon scalar density is recalculated using the new nucleon mass. 
This is repeated at each kp until self consistency is reached when the change in the nucleon 
mass is AM n < 0.1 MeV (the second self-consistent loop). Then, kp is incremented another 
step, and the process continues until one reaches any value of kp desired. A diagram of this 
process is in Fig. 13.21 which clarifies the double self-consistent loop structure. 

3.3 The Numerical Results 

The mass of a free nucleon is computed to be M^{kp = 0) = 1209 MeV. The ~ 30% 
difference is as expected in the model at leading order in Nc- We evaluate the nucleon 
mass Eq. 1)3.15)1 and energy per nucleon Eq. 1)3.11)1 as a function of kp. We choose our free 
parameters to fit E/A — Mtv(0) = B = —15.75 MeV at the minimum. The parameters are 
summarized in Table l3~Tl We use the value g s = 1.27 (corresponding to djv = 41.4 MeV), 




(3.15) 
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Figure 3.2: Diagram of the two-tier self-consistent loop. One inputs a value for fcp and a 
guess for the profile function. O represents the field equation for the profile. represents 
the Dirac equation solved in the finite basis, which takes the profile function as an input, 
and outputs a set of energy eigenvalues E n and wave functions. The E n are then used 
to calculate a the nucleon mass Mjy which is then used to calculate a new nucleon scalar 
density. 
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Table 3.1: Model parameters to fit binding energy B = —15.75 MeV, the location of the min- 
ima, compressibility, and the resulting density of nuclear matter in units of po = 0.178 fm -3 
consistent with the experimental value of 0.167 ± 0.018 fm~ 3 |2U| . 



9s 


9v/^ 


kp 


K 


p/po 


[MeV] 




[fm" 1 ] 


[MeV] 




0.89 


7.22 


1.81 


291.7 


2.26 


1.10 


8.96 


1.51 


312.5 


1.31 


1.27 


10.55 


1.38 


348.5 


1.00 



and vector coupling <?^/47r = 10.55, which gives a Fermi momentum of kp = 1.38 fm -1 in 
nuclear matter consistent with the known value kp = 1.35 ±0.05 fm -1 20 . This value of g s 
is used for the calculation of all of the parton distributions and electromagnetic form factors. 
We plot the binding energy per nucleon using Eq. I)3.11jl in Fig. 13.31 The compressibility 
is K = 348.5 MeV which is above the experimental value K = 210 db 30 MeV, but well 
below the Walecka model |1U8| value of 560 MeV. We have performed the self-consistent 
calculation to obtain the profile functions for zero density, 0.5poi 1-Opo an d l-5po i n Fig. 13.41 
(where po = 0.178 fm -3 is nuclear density, consistent with the experimental value of 0.167± 
0.018 fm" 3 HO]). 

We need to address the valence-only approximation in Eq. l)3.5j) to the self-consistent 
profile solution. Specifically, we are making the approximation that when the change in 
the valence energy level AE V is small, so is AM/v- We are also hypothesizing that the 
negative energy level contributions to the scalar density are approximately constant in 
space, and that the spatially varying part of the scalar density is well approximated by 
the valence level. It is very important to note that these approximations are only used 
in determining the self-consistent profile, not in the wave functions used as input into the 
quark distributions or electromagnetic form factors. We can see the difference from a fully 
self-consistent calculation of the profile function for a free nucleon in Fig. 13.51 The resulting 
profiles are virtually identical. 
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k F [frrT 1 ] 

Figure 3.3: Binding energy per nucleon B = E/A — Mjy as a function of Fermi momentum 
fit to B = —15.75 MeV at the minimum for g s = 0.89 (long dashed), 1.10 (short dashed) 
and 1.27 (solid). The box and shaded region are the experimental uncertainty [20] in the 
binding energy, density and compressibility of nuclear matter. 
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0.5 1 1.5 2 

r [fm] 

Figure 3.4: Profile functions in nuclear matter. The solid line is the profile function for 
1.5/9o; the curves with progressively longer dashes correspond to 1.0po> 0.5po and zero density 
respectively. 
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Figure 3.5: Profile functions for the free nucleon using the valence approximation of Eq. 1)3. 5 J) 
(solid) and a fully self-consistent treatment (dashes). 
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While the vacuum value of the condensate does not vary with the Fermi momentum by 
definition, the effective condensate Eq. (|3.6|) falls ~ 30% at nuclear density, q.v. Eq. ()3.3|) , 
This is consistent with the model independent result [3^ that predicts a value 25-50% below 
the vacuum value. The nucleon mean square radius is plotted as a function of the Fermi 
momentum in Fig. 13.61 The influence of the nuclear medium on the nucleon causes the root 
mean square radius of the baryon density to increase by 3.1% at k F = 1.38 fin -1 . This 
swelling is consistent with a < 6% increase {{r 2 ) a/ {f 2 )o < 1-12 in Fig. I3.6JI as constrained 
by quasi-elastic inclusive electron- nucleus scattering data |75j . 
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Figure 3.6: The nucleon mean square radius as a function of Fermi momentum. The 3.1% 
increase at k,F = 1-38 fm _1 is indicated by the horizontal and vertical lines. Experiments 
|75l ll(J2| limit nucleon root mean square radius swelling to less than 6%, or 1.12 on the 
vertical axis. 
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Chapter 4 
QUARK DISTRIBUTIONS 

One way to probe intermediate length scales and non-perturbative physics is to examine 
the short distance structure of a large object. The prime example is the European Muon 
Collaboration (EMC) effect [Jj where the short distance (~ 5 GeV, or ~ 1CP 2 fm) structure 
of nuclei seen in deep inelastic lepton-nucleus scattering differs from that of a collection of 
free nucleons. This measurement showed that bound nucleons are different than free ones, 
and implied that the medium modifications could be significant for any nuclear observable 

Our primary concern is the depletion of the nuclear structure function F^-(x) in the 
valence quark regime 0.3 < x < 0.7. While the general interpretation is that a valence 
quark in a bound nucleon has less momentum than in a free one, corresponding to some 
increased length scale, the specific mechanism for this has eluded a universally accepted 
explanation for 20 years |50l 1831 H3 I93j . A popular explanation is the so-called 'binding' 
effect which originates from a possible mechanism in which mesons binding the nucleus 
carry momentum. An important consequence is that the mesonic presence would enhance 
the anti-quark content of the nucleus |12l 144] . Such an effect has not been seen in Drell- 
Yan experiments ^2 in which a quark in a proton beam annihilates with an antiquark in a 
nuclear target producing a muon pair. Furthermore, relativistic treatments, including the 
light-cone approach needed to obtain the nucleon structure function, of the binding effect 
with structureless hadrons fail |791 1981 IToT 147] . suggesting that modifications of the internal 
quark structure of the nucleon are required to explain the deep inelastic scattering data. 

Any description of the EMC effect must be consistent with the constraints set by both 
deep inelastic scattering and Drell-Yan data. Thus a successful model must include an- 
tiquarks as well as quarks, and show how the medium modifies both the valence and sea 
quark distributions. Our purpose is to provide a mechanism for that modification within the 



34 



Chiral Quark-Soliton (xQS) model |68 | ll6 | l39 l l26j. This phenomenological model has many 
desirable qualities: the ability to describe a wide class of hadron observables with surprising 
accuracy, the inclusion of antiquarks, positivity of Generalized Parton Distributions, and a 
basis in QCD [22] ■ Here we show how the model describes the EMC effect, and satisfies the 
bounds on nuclear antiquark enhancement provided by Drell-Yan experiments. 

Polarized lepton-nucleus scattering experiments are another important tool in hadronic 
physics. For example, in order to study the spin structure function of the neutron, one must 
use nuclear targets. It is already well known that there are significant differences betweeen 
free and bound nucleons in the unpolarized case (the EMC effect |? ). It is reasonable to 
assume that nuclear effects could appear in polarized quark distributions. Our purpose here 
is to calculate the analogous modification to the nucleon spin structure function gf ,n ^ (x, Q 2 ): 
a 'polarized EMC effect'. 

4.1 Light-cone Correlation Function and Sum Rules 

If we insert a complete set of states into the light-cone correlation function, we obtain 

q[x) = k + [ — e lxk+z ~ (£#(0) 7 + V(O|A;> (4-1) 

J 27T 

= N c k + Y [ — e ixk+z ~ (k\ P ,n}7 j + e- ip " z ~ (p,n\k). (4.2) 

Evaluating Eq. (|4,2j) in the rest frame, with k + = k~ = Mjy[0(r, kp)\ from Eq. (|3.15|) . and 
writing the wave functions and plus and minus components explicitly, we obtain 

q(x) = N C M N £ f ^ V M ^Vt (p)(l + 7 V)e- i(iW)2C V„(p) (4-3) 

We can put this in a form that is ideal for evaluation in the finite KR basis 

q(x) =N c M N Y,{n\{l+ 1 Q 1 Z )5(E n + P i -xM N )\n) (4.4) 
n 

In order to address the baryon and momentum sum rules, we will follow Ref. It is 
convenient to put the light-cone correlator in the form of a functional trace using (in the 
rest frame) 

\n)(n\e- iEnZ ° = J dE5(E - H)e~ iEz ° 
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where H is the Dirac Hamiltonian Eq. (|3.12j) . With this substitution, we can write the 
moments of Eq. (|4.3|) as 

J dxx n - l q{x) = N c M]f n j dETi [5(E - H)(E + p 3 )™" 1 ^ + 7°7 3 )] (4.5) 
For n = 1, we obtain the baryon sum rule 



N c TV 



9{E V - H) - 0(-hW) 



N C B 



since the ^-function operator simply counts the occupied levels in the nucleon, with the 
number of occupied levels in the vacuum subtracted, leaving only the quarks in the valence 
level. It should be noted here that the topological condition for the profile function guar- 
antees that only one energy level comes out of the positive continuum, and hence that the 
winding number of the soliton field is not the baryon number; it is only a necessary, not 
sufficient, condition to obtain unit baryon number. 

For n = 2, along with other terms that vanish by the hedgehog symmetry of the profile 
function, we obtain the momentum sum rule 



N C M^ 1 Tr \o(E v - H)H - 9{-H {0) )H {0) 



NcMZ 1 



E n <0 
M- l M N [<d(r, k F )] 

1 



£< 0) <o 



(4.6) 



The only details we have used are that we have a complete set of basis states, a single 
valence energy level, and that the sum of the energy eigenvalues of the Dirac Hamiltonian 
Eq. 1)3.12(1 yield the nucleon mass Eq. (|3.15|) . Therefore, the calculation presented here is 
independent of the addition of a non-zero vector potential. It is important to note that 
both of these sum rules are automatically satisfied as long as one remains in the same basis 
consistently throughout the calculation, provided Eq. I|4.1|) holds. 



4.2 Unpolarized Distributions 



The isoscalar unpolarized distribution q(x) = u(x) + d(x) is the leading order term in Nq, 
with the isovector unpolarized quark distribution u(x) — d{x) smaller by a factor ~ 1/N C 
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and set to zero. The distributions are calculated using the KR basis at k F = and kp = 
1.38 fm _1 almost exactly as in Ref. jSHl where the quark distribution is given by the matrix 
element Eq. (|4.4j) evaluated in the finite basis 

q(x) = N C M N ^ Yl c *ni c njJ-°(Pi ~ \ E n ~ xM N \) 



E n -xM N Q 

Pi 



(4.7) 

with the Pauli-Villars regulated sum taken over occupied states. The eigenvalues E n are 



2pi 



determined from diagonalizing the Hamiltonian, derived from the Lagrangian (|2.10j) . in the 
KR basis. The vector meson exchange is not explicit in Eq. ()4.7|) ; it is implicit in the energy 
eigenvalues E n . The antiquark distribution is given by q(x) = —q{—x) where the sum is over 
unoccupied states [SHI- We use the exact sea wave functions, and not the approximation 
used in Eq. (J3.5|) . The use of a finite basis causes the distributions to be discontinuous. 
These distributions are smooth functions of x in the limit of infinite momentum cutoff 
and box size, but numerical calculations are made at finite values and leave some residual 
roughness. This is overcome in Ref. jSH] by introducing a smoothing function. We deviate 
from their procedure, and do not smooth the results; instead we find the subsequent one-loop 
perturbative QCD evolution to be sufficient. This can be seen in Fig. 14.11 

These distributions are used as input at the model scale of Q 2 = Mp v ~ 0.34 GeV 2 for 
evolution to Q 2 = 10 GeV 2 in the case of the quark singlet distribution q s (x) = q(x)+q(x) oc 
FJ?(x)/x at leading order in No- We define the generic ratio function to be 

R( x n 2 ) - f®9(x,Q\k F ) 

f®g(x/y,Q 2 ,k F ) = / dyf(y)g(x/y,Q 2 ,k F ). (4.9) 

J X 

In the following, we take the general function g to be the structure function F^ N \ the 
antiquark distribution q, or the spin structure function g^f ,n \ The nucleon momentum 
distribution f(y) following from a light-cone approach, for any mean field theory of nuclear 
matter (see Ref. jTS] and Eq. (|D.40I) in Appendix EI is 

f(y) = ^0(1 + A F - y)9{y — 1 + A F ) [Af- - (1 - y) 2 ] , (4.10) 

where A^ = kp/M^ and = Mtv(0) — 15.75 MeV. The distributions are evolved to 
Q 2 = 10 GeV 2 for use in the EMC ratio Rf 2 {x) and the Drell-Yan ratio R q (x) Eq. (j£K)) . 
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Figure 4.1: The smoothing effect of evolution on the quark singlet distribution. The curves 
are shown at scales log Mp V (highest curve at x = 0.3), and the DIS scale logQ 2 , with 
Q 2 = 10 GeV 2 (lowest curve at x = 0.3). 
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The ratios are plotted in Figs. l4~2l and EOl respectively. While the data shown in Figs. EOl 
and 14. 31 are for large, but finite, nuclei, our calculation reproduces the trend of both sets of 
data. It falls slightly below the SLAC-E139 data [22] due to the higher density of nuclear 
matter. 

In Fig. 14.41 we show the quark singlet distribution for a free and bound nucleon at a 
scale Q 2 = 10 GeV 2 . There is a large depletion in the bound nucleon valence distribution, 
accounting for nearly all of the depletion in Fig. 14.41 since the antiquark contribution is 
virtually unchanged. If one uses only the valence contribution to calculate the EMC ratio 
(|4.8|) . produces too large an effect. This large effect is comparable to that of the QMC 
model impulse approximation calculation or the Guichon model 90 , both of which only 
include valence quarks. This valence effect is mitigated by a small enhancement in q(x), 
consistent with the Drell-Yan data, so that the singlet distribution has only a moderate 
depletion consistent with the EMC effect. 

A simple picture in terms of the uncertainty principle is available. The influence of the 
nuclear medium on the nucleon causes the root mean square radius of the baryon density 
to increase by 3.1%. This corresponds to a decreased momentum, and a depletion of the 
bound structure function relative to the free one. This swelling is consistent with a < 6% 
increase as constrained by quasi-elastic inclusive electron-nucleus scattering data |75| . and 
the recent polarization transfer measurement |102j . 

We ignore the effects of shadowing, which occur when the virtual photon striking the 
nucleus fluctuates into a quark-antiquark pair over a distance ~ 1/2Mnx exceeding the 
inter-nucleon separation. This causes a depletion in the structure function for x < 0.1 
and is relatively well understood [H3 E31 El HHj and so we do not reiterate those results. 
Additionally, we ignore contributions from quantum pion structure functions, which in this 
model propagate through constituent quark loops, and would modify the behavior at small 
x. These loops are suppressed by 0(1/Nc), and are not treated at leading order. 

The present model provides a intuitive, qualitative treatment that maintains consistency 
with all of the free nucleon properties calculated by others |39[ I26j . It gives a reasonable 
description of nuclear saturation properties, reproduces the EMC effect, and satisfies the 
constraints on the nuclear sea obtained from Drell-Yan experiments with only two free 
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Figure 4.2: The EMC ratio at a scale Q 2 = 10 GeV 2 for nuclear matter. The data are for 
Iron (empty boxes) and Gold (filled boxes) from SLAC-E139 |52| . 
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Figure 4.3: The Drell-Yan ratio at a scale Q 2 = 10 GeV 2 for nuclear matter. The data are 
for Iron (empty boxes) and Tungsten (filled boxes) from FNAL-E772 
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parameters: g s and g v . 

The central mechanism to explain the EMC effect is that the nuclear medium provides 
an attractive scalar interaction that modifies the nucleon wave function. This is also the 
dominant mechanism in the QMC model approach to the EMC effect [001 an d also similar to 
the quark derealization approach [0j. The improvements given here are the explicit compu- 
tation of the effects of the medium on the antiquark distributions so that consistency with 
the Drell-Yan data could be verified, and the reduction of the number of input parameters 
and model assumptions. Our extension of the Chiral Quark-Soliton model to nuclear matter 
provides a new, consistent way to calculate possible medium modifications of a variety of 
observables that could be measured in experiments. 

4.3 Polarized Quark Distribution Function 

The first discussion of nuclear effects in the polarized quark distributions is in Ref. [HOI in 
the context of dynamical rescaling. A more recent calculation is in Ref. |28| : it predicts 
dramatic effects for the bound nucleon spin structure function. We have previously shown 
how the model describes nuclear saturation properties, reproduces the EMC effect, and 
satisfies the bounds on unpolarized nuclear antiquark enhancement provided by Drell-Yan 
experiments [001- Therefore, we expect the xQS model to produce a reasonable result for 
the polarized distributions. 

The polarized quark distribution for flavor i is defined by the difference between the 
quark distributions with spin parallel (|) and antiparallel (j) to the nucleon 

A qi (x, Q 2 ) = qj(x, Q 2 ) - q\{x, Q 2 ). (4.11) 

— t — I 

The polarized antiquark distribution is defined analogously using ql , and q\. The isovector 
polarized distribution Aq^ T=l \x) = Au(x) — Ad(x) is the leading order term in Nc, with 
the isoscalar polarized quark distribution Aq( T=0 ' (x) = Au(x) + Ad(x) smaller by a factor 
~ 1 /Nc and set to zero. We will therefore suppress the isospin superscript in the following. 
The distributions are calculated using the KR basis at kp = and kp = 1.38 fm _1 (see 
Refs. ffiH 1100] ) almost exactly as in Ref. [HE] where the quark distribution is given by the 
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matrix element 

Aq(x) 



-\(2T 3 )N C M N V V c^cnj^-eipi - \E n - xM N \) 



i -J 



E n - xM N 
Pi 



(75P • r)ij + 



1 3 (E n -xM N f 



2ft 2 pi 
1 1 (£ B - xM N ) 2 



/ 1- 



(7°r- 7 ), 



(4.12) 



.2 2 pj 

with the Pauli-Villars regulated sum taken over occupied states. The eigenvalues E n are 
determined from diagonalizing the Hamiltonian, derived from the Lagrangian (|2.1()j) . in the 
KR basis. 

The antiquark distribution is given by Aq(x) = Aq(—x) where the sum is over unoccu- 
pied states. The use of a finite basis causes the distributions to be discontinuous. These 
distributions are smooth functions of x in the limit of infinite momentum cutoff and box 
size, but numerical calculations are made at finite values and leave some residual roughness. 
This is overcome in Ref. [20] by introducing a smoothing function. We deviate from their 
procedure, and do not smooth the results; instead we find that performing the one-loop 
perturbative QCD evolution |S7] provides sufficient, but not complete, smoothing. Some 
residual fluctuations due to the finite basis remain visible in our results, and the size of 
these fluctuations serve as a guide to the size of the error introduced by the method. 

These distributions are used as input at the model scale of Q 2 = Mp V ~ 0.34 GeV 2 
for evolution to Q 2 = 10 GeV 2 . The polarized structure function to leading order in Nc is 
given by 



(p> n ) 
9i 



(Pi n ) 
9i 



Aq 



NS 



±-Aq NS + 0(N° C ) 
I (An - Ad + Au - Ad) + O(iVg) 



(4.13) 

(4.14) 
(4.15) 



The nucleon momentum distribution is assumed to be the same as the unpolarized case, 
Eq. (jUmj). 
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This is justified since the Hugenholtz-van Hove theorem follows from the pressure van- 
ishing for a stable nucleus (see Appendix [DJ), which would remain valid for a polarized 
system. Therefore, we would expect to still have a distribution sharply peaked around 
y — 1. The details of the functional form of f(y) have little effect on the ratio beyond the 
peak location. We are already ignoring target polarization, which has the same effect as the 
leading order correction (change in the normalization) to the momentum distribution, with 
the sum on the spin-dependent distributions yielding the unpolarized distributions. Also, 
the spin-orbit force tends to average out in nuclear matter, leading to the same distribution 
for all nucleons. The nuclear matter momentum distribution we use is fairly close to those 
calculated for Lithium isotopes |92j . 

We show the ratio Eq. (|4.8|) in Fig. 14.51 using a 'valence'-like distribution, as well as for 
the full distribution. The latter includes all medium modifications, while the former dis- 
tribution uses the medium modified energy level eigenstate, but the same free nucleon sea 
quark distribution for both the free and bound nucleon. This was done in order to compare 
our results with the model in Ref. |2HJ ? which only has valence quarks at the model scale. 
The single energy level actually has a contribution to the polarized antiquark distribution, 
so it alone cannot be considered a true valence spin structure function. However, this con- 
tribution is small, so we effectively reproduce the result of a valence quark model, especially 
in the region x > 0.3. 

In Fig. 14.51 one can see that there is a large depletion in the polarized 'valence' quark 

(T=l) 

distribution . This produces a large depletion in the isovector axial coupling g A of 
17.8%. This large effect is comparable to that of the calculation in Ref. j^H] which only 
includes valence quarks at the model scale. This valence effect is mitigated by a large 
enhancement in the sea quark contribution, so that the full polarized distribution has only 
a moderate depletion in the region 0.3 < x < 0.7 of the same size as the EMC effect in 
unpolarized nuclear structure functions. There is a large enhancement for x < 0.3 due to 
the sea quarks. This large enhancement is very different from the small effect calculated in 
the unpolarized case [22]) and seen in unpolarized Drell-Yan experiments 2 . This would 
suggest that one might see a significant enhancement in a polarized Drell-Yan experiment, 
even after including shadowing corrections (which we address later). 



45 




Figure 4.5: The ratio Eq. (|4.8|) at scale Q 2 = 10 GeV 2 for nuclear matter. The heavy line 
is the full calculation for nuclear matter. The light line is the effect calulated using only 
medium modifications to the 'valence' energy level as decribed in the text. 
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The larger sensitivity to the lower components of the wave functions, and the smaller 
size of the relevant distribution Aq(x) are the primary sources for the greater sea quark 
enhancement in the polarized case than in the unpolarized case. This model gives a smaller 

(T=l) 

value for the isovector axial coupling g A =1.11 than experiment. This is to be expected 
at leading order in Nc as there are important corrections due to rotational modes of the 
soliton. These corrections should be relatively minor when considering the ratio Eq. (|4.8|) . 

(T=l) 

The axial coupling g A is enhanced by 9.8% in the nuclear medium. This is in accord 
with an earlier finding of a ~ 25% enhancement for gA in a different soliton model by Birse 
|14| . There, the effect is also seen as a competition between enhancement and depletion. 

In order to address the medium modification of the Bjorken sum rule |18[ I19| 



as an integral of the experimentally observed nuclear distribution, one must account for the 
effects of shadowing. This occurs when the virtual photon striking the nucleus fluctuates into 
a quark-antiquark pair over a distance ~ l/2M^x exceeding the inter-nucleon separation. 
This causes a depletion in the structure function for x < 0.1 and is relatively well understood 
|831 R)ll9"3]. Shadowing in the polarized case is expected to be larger than in the unpolarized 
case by roughly a factor of 2 simply from the combinatorics of multiple scattering (see 
e.g. Ref. jSH]). 

The enhancement at x ~ 0.1 — 0.2 in Fig. 14.51 is comparable to that seen by Guzey 
and Strikman pBj; they assume that the combined effects of shadowing, enhancement, and 
target polarization lead to the empirical value of the nuclear Bjorken sum rule for 3 He and 
7 Li. Shadowing effects become large for x < 0.05, but we ignore them as well as target 
polarization; such precision is not necessary for our relatively qualitative analysis. One 
needs ~ 10 times the shadowing observed in the unpolarized case for Lead in order to 
counter the enhancement at x ~ 0.1 — 0.2, and give the same value for the Bjorken sum rule 
(|4.16|) in matter and free space. This assumes that shadowing is the only effect neglected 
at small x in our calculation of the unpolarized quark distribution 99 . 

We also present, in Fig. 14.71 the results for the spin asymmetry 




(4.16) 



A { f\x,Q*) 



E^Ag^Q^ + Aq^Q 2 )] 
E l e 2 [g l (x,g 2 ) + g l ( 3 ;,Q2 ) ] 



(4.17) 
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The nuclear asymmetry A\ is defined by replacing the polarized and unpolarized quark 
distributions, represented generically as q, with / ® q in Eq. (|4.9j) . We find that for the 
free case, the calculation falls slightly below the data due to the smaller value of qa in the 
large Nq limit, and that the size of the medium modification is of the same order as the 
experimental error for the free proton 

The central mechanism to explain the EMC effect is that the nuclear medium provides 
an attractive scalar interaction that modifies the nucleon wave function. We see this again 
in the polarized case. This is also the dominant mechanism in the Nambu-Jona-Lasinio 
(NJL) type quark-diquark model of Cloet et al |28j . and the soliton model of Birse |14| . 

It is worth noting that the model of Ref. [2B] uses proper time infrared and ultraviolet 
cutoffs. These cutoffs are known [HE] in the xQS model to ruin the equivalence of sums 
using occupied and unoccupied states, producing unphysical behavior in the antiquark dis- 
tributions. This lack of equivalence can be traced to a violation of causality, therefore the 
NJL model of Ref. [2Hj, like the MIT bag model used in the QMC model, probably could 
not be used to calculate medium effects on sea quarks. 

The present model provides a intuitive, qualitative treatment that maintains consistency 
with all of the free nucleon properties calculated by others |39| I26j . It provides reasonable 
description of nuclear saturation properties, reproduces the EMC effect, and satisfies the 
constraints on the nuclear sea obtained from Drell-Yan experiments with only two parame- 
ters for the nuclear physics (g s and g v ) fixed by the binding energy and density of nuclear 
matter. Therefore, we expect the results presented here to manifest themselves in future 
experiments with polarized nuclei. 
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Figure 4.7: The asymmetry Af Eq. Q4.17JI at scale Q 2 = 10 GeV 2 . The heavy line is for 
nuclear matter. The dashed line is for the free proton. The data are from SLAC [2] (filled) 
for Q 2 ~ 1 - 40 GeV 2 and HERMES [T] (empty) for Q 2 ~ 1 - 20 GeV 2 The free curve falls 
slightly below the data due to the lower value of gA calculated in the large Nq limit. 
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Chapter 5 

ELECTROMAGNETIC FORM FACTORS 

Recent polarization transfer experiments at TJNAF |102j observed a difference in the 
electromagnetic form factors of a proton bound in a Helium nucleus compared to a free one. 
The form factors encode the structure information of the proton observed in elastic scatter- 
ing. There is extensive work on the medium modifications of electromagnetic properties of 
the nucleon in the literature (for example, see Refs. [571 Hoi ITTU1I73] ). This includes effective 
Lagrangians as well as models that include the quark substructure of hadrons. While in 
principle these effects could be couched in terms of effective field theory operators, it is our 
thesis that such results may be more transparent, physically intuitive or straightforward to 
calculate when viewed as a change in the internal structure of the hadrons. 

We will use the xQS model to calculate the electromagnetic form factors. The overall 
procedure is similar to the Quark- Meson Coupling model (QMC) 73 , which uses the MIT 
bag model for the nucleon. The bag model does not include sea quarks. It is a confining 
model, whereas the xQS model is not. Additionally, the QMC model calculation, when 
coupled with a Relativistic Distorted Wave Impulse Approximation (RDWIA) calculation 
jlU6j or a Relativistic Multiple-Scattering Glauber Approximation (RMSGA) calculation 
|89l I71j , improves the agreement between theory the TJNAF data jlU2j . With our study, 
we hope to reinforce the interpretation of the medium effect in terms of quark degrees of 
freedom, as well as provide an alternate model when the accuracy of the data is improved. 

The procedure is much the same as that in Chapter |I] for the quark distributions, except 
that we are looking at the spatial structure of the nucleon in elastic scattering (after a 
Fourier transform) instead of the momentum distributions of the constituents observed in 
inelastic scattering. In fact, Generalized Parton Distributions (GPD) contain both form 
factors and quark distributions in the same framework. We will begin our description of 
electromagnetic form factors unconventionally starting with GPDs. 
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The leading twist GPD is given by the off-forward (incoming nucleon has a different 
momentum k than the outgoing one k') matrix element analogous to Eq. Q1.8|) 

F(x,tQ 2 ) = k + J g e lxk+z -(k'mOh + i;(z-)\k) (5.1) 

which can be written in terms of two functions H and E 



F{x,i,Q 2 )=u{k') 



^H(x,tQ 2 ) + l ^^E(x,^,Q 2 ) 



u(k) (5.2) 



The functions u(k) and u(k') are the incoming and outgoing nucleon spinors. The other 
variables are defined by 

q = k-k' (5.3) 

Q 2 = -q 2 (5.4) 

^ - q+ (5.5) 



(k + k')+ 

The parton distributions of Chapter |I] are given by the forward limit k = k', so that 

q(x) = H(x,0,0) (5.6) 
Aq(x) = E(x,0,0) (5.7) 

The electromagnetic form factors are integrals over x of functions H and E 

J dxH{x,£,Q 2 ) = Fi(Q 2 ) (5.8) 

J dxE{x,i,Q 2 ) = F 2 {Q 2 ) (5.9) 

where the are the usual Dirac and Pauli form factors 1 . The dependence on £ drops 
out as a special case of polynomiality: the moments of H and E must be polynomials in £. 
This property is rooted in Lorentz invariance, and is preserved by the xQS model |94| . The 
Dirac/Pauli form factors can be rewritten as the Sachs form factors 

G E (Q 2 ) = Fl (Q 2 ) - ^-F 2 (Q 2 ) (5.10) 
G M (Q 2 ) = Fi{Q 2 ) + F 2 {Q 2 ) (5.11) 

lr There are not to be confused with the structure functions F^\x) discussed Chapters Q and 0] 
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These electromagnetic form factors can be evaluated in terms of the quark wave func- 
tions, and are derived in Ref. (27]. The formulae are reproduced here, with a Pauli-Villars 
regulator, for convenience. To leading order in Nc, we have only the isoscalar electric and 
isovector magnetic form factors (G x ~ 0,1 = G x ± G x ) 

G T E =°{q 2 ) = ^fdre^A £ ^ (r)^(r) - £ 4°)t(r)^°)(r)l (5.12a) 
GlT\ q 2 ) = ^L £ ^ fdre^l £ V^W/A^r)} 

° I? M [ E„<E V J 

-PI/. (5.12b) 

In the nuclear medium, Eqs. 1)5.12(1 acquire a dependence on the Fermi momentum through 
the wave functions as described in Chapter El We use Eqs. ()5.12aJ) and (|5.12b|) to calculate 
the form factors, which we present in Figs. 15.11 and 15.21 We also present the results in 
terms of the ratios 

G T E :»j l (Q\k F ) _ G X (Q 2 ) 

g t e :$Hq 2 ,o) -G X m' [b - 6) 

where X is E(T = 0) or M(T = 1), and the double ratio 

G E (Q 2 )/G* M (Q 2 ) 

Ge{Q 2 )/Gm(Q 2 ) ' K ' ' 

These ratios are plotted in Figs. 15.31 and 15.41 for 0.5po, 1-Opo and 1.5/?o- 

The electric form factor is dominated by the valence contribution and shows a dramatic 
effect, while the magnetic form factor has equally important contributions from the valence 
and the sea. The latter shows almost no change in nuclear matter; it shows only a 1.3% 
enhancement of the magnetic moment at nuclear density, and a 2.3% enhancement at 1.5 
times nuclear density. These are consistent with the constraints of a < 2% increase set 
by Ref. |45| . In contrast, the QMC model |73j predicts an enhancement of roughly 5% for 
Oxygen and 10% for Lead. The influence of the nuclear medium on the nucleon also causes 
the root mean square radius of the baryon density to increase by 3.1% in the xQS model. 
This swelling is consistent with a < 6% increase as constrained by quasi-elastic inclusive 
electron- nucleus scattering data 75 . The effect in the electric form factor calculated here 
is comparable to that of the QMC model; the main difference from that calculation lies in 
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Figure 5.1: The isoscalar electric form factor at nuclear density (solid) and at zero density 
(dashes). 




Figure 5.2: The isovector magnetic form factor at nuclear density (solid) and at zero density 
(dashes). 
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Figure 5.3: The electric (lower three curves) and magnetic (upper three curves) form factor 
ratios in Eq. (|5,13|) for 0.5po (long dashes), I.O/jq (solid) and 1.5po (short dashes). 
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0.5 1 1.5 2 2.5 

Q 2 [GeV 2 ] 

Figure 5.4: The double ratio Eq. Q5.14|) of the electric to magnetic form factors in nuclear 
matter and in the vacuum from the xQS model (heavy) and the QMC model [73] (light). 
Three densities are shown: 0.5po (long dashes), 1.0/?o (solid) and 1.5po (short dashes). 
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the lack of enhancement in the magnetic form factor, specifically the practically unchanged 
value of the magnetic moment. 

While both form factors use the same wave functions, the isovector magnetic form factor 
includes an extra weighting by a factor of the angular momentum of the state (relative to the 
electric form factor) due to the -y k in Eq. (|5.12b|) . This extra factor is not only responsible 
for making the regularization of Eq. ()5.12b|) necessary, but for making the sea contribution 
as important as the valence. In the xQS model, the orbital angular momentum carried by 
the sea is comparable to the orbital angular momentum carried by the valence quarks |107j 
(the sum of which make up about 60% of the total angular momentum of the nucleon state, 
with the remainder belonging to the intrinsic spin of the constituent quarks). 

Conversely, the isoscalar electric form factor (which is finite, after the vacuum subtrac- 
tion) does not have as large of a contribution from the sea. The valence level is the most 
important piece, even at Q 2 > 0, since the Q 2 dependence in the form factors arises from 
the wave functions |27| . The negative Dirac continuum wave functions largely cancel in the 
vacuum subtraction in Eq. (j5,12aj) . 

The magnetic form factors are sensitive to the tail of the quark wave functions, and 
the mere existence of a tail is due to the lack of confinement. This is one reason for the 
discrepancy between the current results and the QMC model [73], but the primary source 
is due to the resistance to change of the sea. The former accounts for only a few percent 
of the difference; it is the latter that is our most important result. We see that the role of 
antiquarks is again prevalent as in our previous work |99j . 

The double ratio obtained in Fig. 15.41 has the same trend as the QMC model [73], but 
differs in the details. Since we obtain a similar double ratio, we expect to have similar 
results if we compare these results with the polarization transfer data |102| . This requires 
one to take the final state and relativistic effects into account through the use of the RDWIA 
|1U6| or the RMSGA which accounts for a few percent of the discrepancy between the 
results for bound and free protons. A RMSGA calculation for the Helium reaction studied 
in Ref. |102| has been done with these xQS model results [ZD]) and it delivers remarkably 
similar results to the same calculation done with the QMC model |71j . The xQS model 
predicts a smaller deviation than the QMC model from a Relativistic Plane Wave Impulse 
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Approximation (RPWIA) calculation, which is taken as a baseline in Ref. |102| . While it 
slightly worsens the agreement with the data at Q 2 < 1, the differences are of the same 
order of magnitude as the current experimental error, and both models under predict the 
observed deviation from a RPWIA calculation. At higher Q 2 , the two models produce 
nearly identical results for Helium. 

We ignore potentially important corrections that follow from integrating out the rota- 
tional and translational zero modes of the soliton that are suppressed by \/Nc- These 
corrections break the N — A degeneracy, and improve the agreement of the vacuum form 
factors with experiment |27j . More relevant to the calculation presented here, the rota- 
tional corrections do not affect the Q 2 dependence, but instead affect the normalization 
of the form factors [27j. However, there is no reason at that level to continue to ignore 
quantum fluctuations of the the pion field (quark loops, also suppressed by 1/Nc) or center 
of mass corrections, and treat the soliton as a purely self-consistent, static mean field. We 
will save this difficult problem for the future. 

We have calculated the electric and magnetic form factors at leading order in Nq at 
nuclear density using the xQS model. Our results help validate the apparent success of 
the QMC model in describing the polarization transfer experiment HQ2J [73|, and provide 
a counterpoint consistent with constraints on the nucleon radius and magnetic moment in 
nuclei to be distinguished when finer resolution becomes available in the data. In fact, the 
difference between the xQS model double ratio and the QMC model [ZS] is roughly the size 
as the current experimental error. Specifically, data on the bound nucleon magnetic form 
factor at low Q 2 , particularly the magnetic moment, could serve to determine the role of 
sea quarks in nuclei. 
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Chapter 6 
SUMMARY 

The Chiral Quark-Soliton model provides a intuitive, qualitative treatment that main- 
tains consistency with all of the free nucleon properties calculated by others |391 126j . It 
gives a reasonable description of nuclear saturation properties, reproduces the EMC effect, 
and satisfies the constraints on the nuclear sea obtained from Drell-Yan experiments with 
only two free parameters: g s and g v . 

The central mechanism behind these calculated effects is that the nuclear medium pro- 
vides an attractive scalar interaction that modifies the nucleon wave function. This is the 
same mechanism involved in other calculations |90 | l9|l28 l ll4j. The major difference between 
this calculation and others is the inclusion of antiquarks at the model scale, and we see that 
they play an important role in the modifications, maintaining not only consistency with 
constraints such as those imposed by Drell-Yan experiments [2] and the theoretical and 
experimental understanding of magnetic moments |45| . but also positivity and baryon and 
momentum sum rules. The latter is deeply connected with the inclusion of antiquarks and 
a regulator which does not interfere with the completeness of the set of quark states or the 
equivalence of summing over occupied and unoccupied states in the spectrum. 

Dynamical rescaling, a moderately successful attempt at simultaneous phenomenological 
description of the EMC effect and related Drell-Yan experiments, automatically maintains 
the baryon and momentum sum rules due to the structure of the renormalization group 
equations. A change in the nucleon size in nuclei was believed to be the scale driving the 
effect. These are all properties of the xQS model, with the nucleon size swelling 3.1% in 
the nuclear medium. Although there is no special scale (one could point to the relative 
sizes of the nuclear scalar density and chiral condensate), we have a mechanism in the 
scalar interaction of the quarks in a nucleon with the medium. It may be possible that the 
spectrum of quark states in a bound nucleon is the exactly the same as in a free nucleon at 
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a smaller model scale Mpv, the Pauli-Villars mass, and hence one would have a 'derivation' 
of dynamical rescaling. It also may just be a coincidence that any model that maintains 
the momentum and baryon sum rules with some swelling of the nucleon describes the EMC 
effect and Drell-Yan data within the accuracy of current experiments. 

The Wilson coefficients c in Eq. ()1.9|) in the OPE are dependent on the starting renor- 
malization scale. Nuclear interactions are effective QCD interactions which renormalize a 
free nucleon into a bound one. In the handbag diagram Fig. 11.21 all three quarks, not just 
the spectator quarks, in the intermediate state are renormalized by QCD from the same 
starting scale, and hence all feel the effects of nuclear interactions. The picture presented in 
Chapter 0] is consistent with this viewpoint, while other approaches |28| I90j treat the struck 
quark as if it had the same normalization point as one in a free nucleon. 

We have also calculated the electric and magnetic form factors at nuclear density. Our 
results help validate the apparent success of the QMC model in describing the polarization 
transfer experiment [1021 FZBL and provide a counterpoint to be distinguished when finer 
resolution becomes available in the data. In fact, the difference between the xQS model 
double ratio and the QMC model [23] is roughly the size as the current experimental error. 

These successful descriptions of the medium modification of quark distributions and form 
factors give weight to the model and point to future experiments. For example, data on the 
bound nucleon magnetic form factor at low Q 2 , particularly the magnetic moment, where 
the difference between the QMC and xQS models is the largest, may be able to differentiate 
between the two models. We also expect the results for the polarized quark distribution to 
manifest themselves in future experiments with polarized nuclei. These experiments could 
serve to help determine the role of both quarks and antiquarks in nuclei. 
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Appendix A 

THE CONNECTION TO CHIRAL PERTURBATION THEORY AND 

NONPERTURBATIVE QCD 

In this Appendix, we perform a low-momentum expansion of the xQS model Lagrangian 
1)2.10)) with a non-zero current quark mass, and compare it to the Lagrangian of Gasser and 
Leutwyler. This has been done before in Ref. [S3]) f° r example; we present a slightly different 
approach here for completeness. We also fit quenched lattice data for the dynamical quark 
mass to a sum of propagators (as appearing in Pauli-Villars regulation of divergences). 

A.l The Chiral Lagrangian 

We will first derive Eq. 1)2. 17|) with current quark mass mo 7^ 0. Then we will proceed to 
higher order in the derivative expansion. 

A. 1.1 The 0(p 2 ) Lagrangian 

The constituent quark model Lagrangian in Euclidean space with (anti)quark fields ip\ip 
that we use is 

£ = ^{ip + iMlP 6 + im )i> - ip\ip + iM)ip (A.l) 

where 

[775 = U^—^- + U^^~y^ (A.2) 

The constituent quark mass, M, is a function of momentum, falling to zero as the momentum 
goes to infinity. The second term represents the vacuum subtraction. One can integrate 
over the constituent quark fields in the Euclidean partition function to obtain the effective 
action 

S[U] = -N c \ logdet(i^ + iMU^ 5 + im ) - logdet(i^ + iM)\ (A.3) 
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We wish to look at the real part of this effective action. 

Nr r 1 + 1 

Re S[U] = — y logdet |l -^-^-(MpU' f5 - m M[U + C/ f ])J (A.4) 

To simplify our calculation we introduce the 'regulated' derivative and scalar source opera- 
tors 

M 

= (A.5a 

-S 2 + M 2 

M 

X = 2m — (A.5b) 

-d 2 + M 2 y ' 

Note that V^U is leading 0(p) and % ls leading 0{p 2 ) from the current quark mass. After 

rewriting our effective action in terms of this new notation, and using log det = Tr log where 

Tr means a functional as well as Dirac and flavor trace, we obtain 

N r r 1 + +1 

*&S[U] = — ^-Tr log [l - flT* - -(x j U + tfx)\ (A.6) 

Then one can expand in derivatives of the field V^t7 and current quark mass, and obtain 
the momentum expansion of chiral perturbation theory. One essentially expands 

log(l - e) = -e - X -e 2 - ^e 3 + ■ ■ ■ 

To order 0(p 2 ), after taking the Dirac trace, we obtain the chiral Lagrangian in leading 
order in Nc 

C (N ,2) = NctrV^UV^ + N c tv( X ] U + U ] x) (A.7) 

Where the lowercase trace is over flavor indices, and functional traces over the operators 
defined in Eq. (|A.5|1 . 

A. 1.2 The C(p 4 ) Lagrangian 

To order C(p 4 ) and leading order in Nc, we obtain the operators 

c {NcA) = ^ti(y»uv u utv v uv^ + v»uv^v v uv u ui 

+Nctr(v fl UV tl U\x ] U + 

N c 



+^tr((x t C7 + U^x?) (A. 
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Table A.l: Constants appearing in the integrals for Li 
i= 1 2 3 5 8 



\ . I I 1 i I 

A « 4 2 1 1 2 

a, 4 4 4 3 2 



We will use the identity, for Nf = 3 

tr(ABAB) = -2tr(A 2 B 2 ) + hi(A 2 )tr(B 2 ) + tr(AB) 2 (A.9) 

to put this Lagrangian in the form of Gasser and Leutwyler |49| 

C^c^ = ^(trV^UV^) 2 + ^tv(V^UV u U t )tv(V ll UV I/ U t ) 
-N c tv(V^UV^V u UV u U^) 
+N c tr(v fl UV^(x ] U + rfx) 



+^v(^Ux 1 U + U^xU ] x) (A.10) 

This means that to leading order in Nq, the constituent quark Lagrangian Eq. (|A.1|) gives 
the Gasser and Leutwyler coefficients L\, L2, L3, L5, and L% 7^ and L4 = Lq = L7 = 0. 
When one takes the functional trace over the operator, one gets 

L^^ipLi *w r ( A . u) 

where Aj and a« are given in Table lA~Tl 
A. 2 The Quark Propagator 

The (quenched) lattice data [SB] gives the quark propagator in terms of functions Z(k), the 
wave function renormalization, and M(k), the constituent quark mass. We fit this data 
with a sum of propagators as the parameterization (much like what appears in Pauli-Villars 
regularization) 

Z{k) fzA , 1 . (A.12) 

^— ' 11.2 ik A- 111: -A- i.m.r. 



iji + iM(k) ^ /ij ijt + i/Uj + im 
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Regularizing the divergent integrals appearing in the momentum expansion (specifically 
the quadratically divergent chiral condensate and the logarithmically divergent pion decay 
constant in C^), along with the restriction that Z(k — > oo) — > 1, requires 

n 

Y, z i— = 1 ( A - l3a ) 

i=0 ^ l 
n 

= (A. 13b) 

i=0 

n 

^Z^i = (A. 13c) 

i=o 

We have taken n = 3, so we have 5 parameters to fit the lattice data. A fit of the momentum 
dependent quark mass M(p) with Eq. (jA.12|) for n = 3 to the quenched lattice data for 
m q ~ 0.14 GeV are shown in Fig. IA.ll The same parameters do a good job of reproducing 
the lattice data for several other quark masses, so the extrapolation to the chiral limit shown 
in Fig. 12.31 is expected to allow a good qualitative comparison with the instanton model and 
the single Pauli-Villars subtraction used in the self-consistent calculations. 
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Figure A.l: Fits of Eq. (jA.12|) for n = 3 to the quenched lattice data |86j for the dynamical 
quark mass M(p) at current quark masses m q a = 0.012, 0.018, 0.024, 0.036, 0.048, 0.072, 
0.108, and 0.144, where a is the lattice spacing. 
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Figure A. 2: Fits of Eq. ()A.12|) for n = 3 to the quenched lattice data for the wavefunction 
renormalization function Z(p) at various quark masses. The data falls slightly below the 
fits. This can be attributed to an effect of the finite volume of the lattice |86j . 
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Appendix B 
THE KAHANA-RIPKA BASIS 

The Kahana-Ripka basis is defined as the following set of functions inside a spherical 
box with radius L: 

<PnGjiMst{r) = N n) iji(p n r)E% jlst (n) (B.l) 
where ji(k n r) is a spherical Bessel function. The momenta p n are defined by 

iciPnL) = 

where p n < A and A is the basis momentum cut-off, and G is the grand spin, defined below 
in Eq. (|B.4|I . The normalization N n i is defined by a property of the zeros of spherical Bessel 
functions. If a and a' are two zeros of ji(x) 

jl(a) = ji(a') = 

then 

dyy 2 ji(ay)ji(a'y) = / dy y 2 ji±i{ay)ji ±1 {a'y) (B.2) 

Jo 

= Sa,a>l[3i±i(a)] 2 (B.3) 

The quantum number n is a good quantum number for the free Dirac Hamiltonian Eq. (|3.13|) . 
The Dirac Hamiltonian in the xQS model Eq. (|3.12|) derived from the Lagrangian Eq. H2.1U|) 
commutes with the grand spin G, the sum of the spin, angular momentum and isospin 
operators 

G = J + T = L + S + T, (B.4) 

therefore we construct eigenstates H of G. The functions H are spinor-isospinor-spherical 
harmonics, defined by the sum of Clebsch-Gordan coefficients Cj^ ni j 2m2 , spherical harmon- 
ics Y and (iso)spinors % 

33 *3 h s 3 



L 
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The functions Eq. (jB.l|) are organized into eight (four in the case of G = 0) basis functions: 



four with parity P = ( — 1)' = (— 1) G and four with P = (—l) 1 = (— 1) G+1 . The eight can 
also be separated into two groups of four with I = G + 1,G (x2) and G — 1. There are 
only four functions for G = since I = G — 1 and j = G — 1/2 cannot be satisfied. The H 
are eight element vectors, corresponding to the four of the Dirac spinor and the two of the 
isospinor. For an explicit example, the basis states with G p = + are 



with N n j = V2/ji +1 ( Pn L). They satisfy I = G = j - 1/2 and I = G + 1 = j + 1/2 
respectively. 

An explicit example of the momenta for a Kahana-Ripka basis with cutoff A = 600 MeV 
and box size L = 5 fm is in Table iBTTl These represent a finite basis of 1576 functions. The 
largest basis used to calculate the numerical results presented here has cutoff A = 5000 MeV 
and box size L = 6 fm, for which there are ~ 1.7 x 10 6 functions. Symmetry considerations 
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Table B.l: The momenta for each grand spin G for the 1576 Kahana-Ripka basis functions, 
with cut-off and box size A = 600 MeV and L = 5 fm, respectively. The number of basis 
functions for each G for A = 600 MeV (I), and A = 5000 MeV (II) are shown in the far 
right two columns. 



G 


Pn=l 


P2 


P3 


PA 


Total states 






[MeV] 


[MeV] 


[MeV] 


[MeV] 


(I) 


(II) 





123.984 


247.968 


371.953 


495.937 


16 


192 


1 


177.334 


304.880 


430.336 


555.128 


96 


1128 


2 


227.457 


358.938 


486.330 




120 


1880 


3 


275.782 


411.116 


540.598 




168 


2576 


4 


322.928 


461.939 


593.546 




216 


3312 


5 


369.321 


511.729 






176 


3960 


6 


414.893 


560.700 






208 


4680 


7 


460.049 








120 


5280 


8 


504.793 








136 


5984 


9 


549.194 








152 


6536 


10 


593.302 








168 


7224 












1576 


42752 



immediately reduce this number to ~ 20, 000 depending on the specific matrix elements 
calculated. 
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Appendix C 
DGLAP EVOLUTION 

Quantum Chromodynamics tells us how to calculate parton distributions at another scale 
Q 2 if they are given at a scale Qq through the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) evolution equations |421 1531 H| . In Chapter we make use of QCD evolution to 
not only obtain the (un)polarized distributions at Q 2 = 10 GeV 2 where DIS takes place, 
but to smooth out the effect of the finite KR basis. First, we must write down the polarized 
distribution in its nonsinglet form, and the unpolarized distribution in its singlet form; these 
are more convenient for the application of the DGLAP equations. 

The unpolarized singlet distribution is given in terms of the distributions of up and down 
quarks as 

q S (x) = + fc(x) (C.l) 

i 

= u(x) + d(x) +u(x) + d(x). (C.2) 

The polarized nonsinglet distribution is 

Aq NS (x) = X ^rffi (A qi (x) + Aq t (x)) (C.3) 

= - (Au(x) - Ad(x) + Au(x) - Ad(x)) + O(N^). (C.4) 

Note that Eqs. (|G2|) and (|C.4|) are normalized to the baryon number and axial coupling, 
respectively, which are both 0(Nc)- The combinations u — d and An + Ad are normalized 
to the isospin and spin, respectively; these are 0(N®). The gluon distribution is denoted 
g(x,Q 2 ). 

If we define 

2 a s (Q 2 ) 
l ~ Po g a s (Q 2 ) 
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and 



- 1 dy 

p®q(x) = I — p{x/y)q{y) 

y 



then the DGLAP equations are given by 

d . , 

where 



P M ® Aq NS (x,t) 
P<g>q(x,t) 



(C.5) 
(C.6) 



( •; 









2iV f P, 



99 



f r qg 



Pc 
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Pc 
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The splitting functions jftj to leading order in a s are 



(1^)7 + 2 5(1 " x) 



T fl [x 2 + (1 - x 2 )] 



C F 
2C G 



l + (l-x) 



X 
X 



(l-x). 



1 - X . 
H h x(l — x) 



11 1 N f T R 

12 3 C G 



5(1 -x) 



(C.7) 
(C.8) 
(C.9) 

(CIO) 



where the plus-distribution is defined by its integral with an arbitrary function /(x) 



(l-x). 



(l-x) 



The DGLAP evolution equations are solved here with the brute force method following 



Refs. 



with Nt = 200 points in the scale variable t, Atj = tj+\ — tj, and N, 



1000 



points in x, Axk = x^ — x^-i- We then have, for the nonsinglet polarized distribution 

N x . 

Aq NS ( Xl ,t j+l ) = Aq NS (x i ,t j ) + Atj —Pqqixi/x^q^ix^tj). 

, Xf~ 
k=i 



and similarly for the singlet unpolarized distribution. 
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Appendix D 

THE FAILURE OF RELATIVISTIC NUCLEAR MEAN FIELD 

THEORY 

We have shown |79| I98| that the use of the relativistic mean field approximation, and 
the assumption that the structure of the nucleon is not modified by effects of the medium, 
to describe infinite nuclear matter leads to no appreciable binding effect. The failure was 
encapsulated in terms of the Hugenholtz-van Hove theorem |H2] which states that the average 
nuclear binding energy per nucleon is equal to the binding energy of a nucleon at the top of 
the Fermi sea. The light front version of this theorem is obtained from the requirement that, 
in the nuclear rest frame, the expectation values of the total plus and minus momentum 
are equal. The original version of the theorem was obtained in a non-relativistic theory in 
which nucleons are the only degrees of freedom. Here, the mesons are important and the 
theory is relativistic, but the theorem still holds. This theorem can be shown to restrict [7S] 
the plus momentum carried by nucleons to be the mass of the nucleus, which in turn implies 
that the probability for a nucleon to have a plus momentum k + is narrowly peaked about 
k + = Ma/ A = Mjy. Thus the only binding effect arises from the average binding energy, 
which is much smaller than the average separation energy. Therefore dynamics beyond the 
relativistic mean field approximation must be invoked to explain the EMC effect. 

We briefly outline our procedure. In Sections ID . II and ID . 21 we present the covariant deep 
inelastic scattering formalism of Ref. [HE] and derive its representation in terms of nucleon 
single particle wave functions. The plus momentum distribution follows from this represen- 
tation in Section fD . 31 where we also derive a relativistic version of the Hugenholtz-van Hove 
theorem. Then we present the results of analytic and numerical calculations in Section ID. 41 
the latter giving an A dependence of the ratio function contrary to experimental results. 
This demonstrates that the use of the relativistic mean field approximation, combined with 
the assumption that the nuclear medium does not modify the structure of the nucleon, 
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Figure D.l: Feynman diagram for deep inelastic lepton- nucleus scattering. A nucleus of 
momentum P is struck by a virtual photon of momentum q. We label nucleon momentum 
k, and quark momentum p. 



cannot describe the EMC effect. The reasons for the subtle differences between the results 
for finite nuclei and nuclear matter are detailed in Section fD. 51 

D.l Nucleon Green's Function for Finite Nuclei 

The use of a manifestly covariant formulation to derive the expression for the structure 
function leads to a convolution formula Eqs. (|1.11|) and 1)1.12(1 |66j . The function x(k,P) is 
(proportional to) the connected part of the nuclear expectation value of the nucleon Green's 
function, and the trace is over the Dirac and isospin indices. We begin with the covariant 
plus momentum distribution function. Eq. (|1.12j) where we identify 

X A (k,P) = -i j d 4 x J d i ye- ik < x - y) G c (x,y), (D.l) 

where G (x,y) is the connected part of the nucleon Green's function: 

iG(x,y) = {P\T + {^{x)^\y)}\P). (D.2) 

This result is directly determined from the Feynman diagram in Fig. lD.ll following Ref. |66| . 
but with x having a different normalization. So far this is independent of the particular 
relativistic mean field model, but for concreteness we use a Quantum Hadrodynamics (QHD) 
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+ • • • 



Figure D.2: The function x in a mean field theory with scalar (er) and vector (u) mesons. 



Lagrangian I95j . specifically QHD-I as in Ref. [2^) where the nucleon fields, that 
appear in Eq. (|D.2|) are those appearing in the Lagrangian. Light front quantization requires 
that the plus component of all vector potential fields vanishes, and this is obtained by using 
the Soper-Yan transformation llllj 

V(x) = e"^ A(:r) *(x), d + A{x) = V + {x). (D.3) 

to define the nucleon field operator for various models |79j . This transformation allows 
the use of the eigenmode expansion for the fields which have been obtained previously in 
Ref. [21] 



= ^2[a a ^ Q (x)+biip c 

a 

= Yl [a a ^a(x)e- i P- x+ / 2 + bii Pa {x)e i P- x+ / 2 ] , (D.4) 

a 

where a a and 6« are (anti-) nucleon annihilation operators and we define z = —x~/2 with 
d + = 2d- = —d z and (a;_|_,z) = x which allows us to treat the minus and perpendicular 
coordinates on equal footing. The ^ a and (p a are coordinate space 4-component spinor 
solutions to the light front Dirac equation with eigenvalues p~/2 = M/y — s a . To simplify 
the analysis we will temporarily ignore electromagnetic effects, but we will include them 
in the final numerical results. The light front mode equations in QHD-I are obtained 
by minimizing the P~ operator (light front Hamiltonian) with the constraint (2^ that 
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P + = P . The result is 
-id z tB~(x) 

with 



[a± ■ (p± - g v V±) + /3(M N + g a </>)] V+(x) 

[-id, + 2g v V~] 9+(x) 

+ [a x • (PL - 9vV L ) + 0{M N + g s <t>)] *~(sb) 



1 



z 



(D.5) 
(D.6) 

(D.7) 
(D.8) 



Using standard manipulations and defining as the energy of the highest occupied 
state, we find the Green's function to be 



G(x,y) = ^ Q (z)* a (y)e-^[A(*)-A(y)] 

a 

+2m5(k~ - p~)6(ef - e a ) 
= G D (x,y) + G c (x,y), 



dk~ 



-ik-(x+-y+)/2 



k — p a + it 



(D. 



where the superscripts D and C represent the disconnected and connected parts of the 
nucleon Green's function, respectively. The connected part is relevant to deep inelastic 
scattering and is given by 



G c (x,y) = iJ2 *a(jc)* a (y)e- <ft ' 



-ig v [A(x)-A(y)] -ip a (x+-y+)/2 



(D.10) 



where the sum is over occupied levels a in the Fermi sea F. We now substitute Eq. (ITXTol) 
into Eq. (|D.1|) . first defining (k±, k + ) = k where k • x = k± ■ x± + k + z = k± • x± — k + x~ /2, 
dx = d 2 x±dz, dk = d 2 kj_dk + and 



We find 



K( k ) = j dxe-^e-^^^aix). 



X A (k,P) = (2vr) 2 ^ ^(k)t^k)6(k- - p-). 



(D.ll) 



(D.12) 



The motivation for the 'double-prime' notation is the subject of the next section. 
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D.2 Wave function Subtleties 

It would be useful to express P) in terms of solutions of the ordinary Dirac equation, 

because one may use a standard computer program jHUl- If we use Eq. (|D.11|) we find that 
these 'double-primed' fields satisfy another version of the mode equations Eq. (|D.5|) and 
Eq. (|D.6|) following from an application of the Soper-Yan transformation Eq. (|D.3|) . and are 
given by 

[~id z - g v V + ] = [a ± ■ (p ± - g v V ± ) + 0{M N + g s ^)} ^+ (D.13) 

[id z + g v V + - 2g v V~ + p~] VP^+ = [a ± ■ (p ± - g v V ± ) + p(M N + g a </,)] <".(D.14) 

If one multiplies Eq. (|D.13|) by 7+ and Eq. (jD.14|) by 7" and adds the two equations, using 
V + = V" = V" = V°, one obtains 

H 3 (id z + p-/2) + 7 °(p-/2 - g v V°)] = [ 7 ± • P± + M N + grf] V^x) (D.15) 

Which looks nearly like the ordinary Dirac equation which the exception of the r y 3 p~/2 
term. To remove this term, set 

* a (x) = e- ip « z/2 ^{x) (D.16) 

and substitute into Eq. (|D.15|) . so that 

7 °(p-/2 - g v V°)V a (x) = [ 1 .p + M N + g s <t>) V a (x) (D.17) 

which looks superficially like the ordinary Dirac equation for the VP' fields which appear 
in the Lagrangian. There is a subtle difference; since the light front energy is given by 
k 2 = (k 2 ± +M^)/k + , Eq. (lP~l"5|l has support for k + > 0. This means Eq. (|D.17|) has support 
for A; 3 > — p^./2 which in turn implies that it can only be considered the ordinary Dirac 
equation with a momentum cutoff at p~/2 ~ M^. This restriction is nearly superfluous 
since the probability that a nucleon is carrying | k 3 | > Mn is suppressed by a factor of order 
with R being the nuclear radius. This allows us to effectively identify VP Q ~ VP' 
which gives the approximate relationships between the (equal time) Lagrangian fields VP', 
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the Soper-Yan transformed light front fields VP and the untransformed light front fields ty" 
that appear in Eq. CEE2J). We have 

V' a (x) ~ e~ ip " z / 2 ^{x) (D.18) 

= e -ipaZ/2 e -ig v A( X )y a ^ p_ 19 ) 

Eq. (|D,19|) is the approximate relationship between the VP and VP' fields in Ref. [U- We now 
are ready to derive a representation of Eq. 1)1.12(1 in terms of these nucleon wave functions. 

D.3 Derivation of the Plus Momentum Distribution 

D.3.1 General Case 

In Ref. |2j, it was determined that a plus momentum distribution in QHD-I is given by 

/(*+) = 2 £ [ d 2 x L \*t(x±,k + )\ 2 . (D.20) 

a£F J 

This distribution peaks at k + /Mjv = y — 0.8 for 16 0, (with smaller values for heavier 
nuclei) but is not the distribution obtained from the covariant formalism of Section ID. II 
The connection between this f(y) and the covariant f(y) was made in Ref. |79j : it was 
determined that, in the limit of infinite nuclear matter, the relationship between f{y) and 
f{y) is simply a shift in the argument by the vector meson potential: 

f(y) = f{y + g v V + /M N ). (D.21) 



This shift arises from the use of the Soper-Yan transformation Eq. (|D.3|) where the fields 
are those appearing in the Lagrangian and are used to determine f(y), whereas the ^ fields 
are used to determine f{y). In finite nuclei, this relationship is somewhat more complicated 
since the vector meson potential is no longer a constant over all space. We start with 



Eq. (|D.12|) . and see that 

Tr 1 + X A {k,P) = (27r) 2 Y J ^[l + K(k)K(k)]5(k--p- 



87r 2 J2K + (k)\ 2 S(k--p-). 
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Substituting into Eq. (|1.12|) we obtain 

/GO = 7A2 E f dk5{y-k+/M N )\K + (k)\ 2 . (D.22) 

Use of Parseval's identity and integrating over k+ gives us our main result: 

f(y) = 2M N Y, [ d 2 x ± \^ + {x ± ,M N y)\\ (D.23) 

so the plus momentum distribution is related to Fourier transform of the W wave functions. 
One can see the similarity to Eq. ()D.20|) ; the difference lies entirely in Eq. (|D.11|) . It should 
be emphasized that this result does not depend on the approximation in Section ID. 21 

We shall use f(y) to compute the nuclear structure function F^{x) in Section fD. 41 but 
first we derive a version of the Hugenholtz-van Hove theorem valid for finite nuclei. To do 



that, multiply Eq. (jD,22|) by y and integrate 



{y) = I dyyf(y) 

2 ^ r „ k+ 



£ j dkt-\K + {k)\\ (D.24) 



(2vr) 2 

Now remove the plus projections and re-express ^f" and its complex conjugate in coordinate 
spaces x and x' . One can then integrate over k yielding a delta function 5 (x — x') which 
allows integration over x' yielding 



We wish to look at the fields in order to understand our result in the context of Ref . [2j| , 
so we need to perform the Soper-Yan transformation Eq. (|D.3|) and use 

<f> = T^E f dx* a (x) 1 + [id + +g v V + (x)}V a (x) 

If we explicitly put in the the nuclear state vectors, we can perform the sum on a by 
inserting creation and annihilation operators; we can add the time dependence for free since 
it is unaffected by d + and cancels with both fermion fields, and the vector potential is static. 
We have effectively undone the substitution Eq. (|D.4|) and now have an expectation value 
of an operator 

(y) = ^ J dx(P\^ + [id + + g v V + ]*\P) (D.25) 
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The second term of Eq. (jD.25|) was essentially included by Birse ^2] via a kinematic argu- 
ment; here it follows from a fully covariant light front treatment. Using the vector meson 
field equation in QHD-I 

+ m l V + = g v * 7 +q,, 

integrating by parts, and anti-symmetrizing one can re-express the second term of Eq. ()D.25|) 

(y) = — [ dx(P\ty-f + id + ty + m 2 v V + V + + V +fl V U + \P) 
Ma J * 

= J dx(P\T ++ - d + (f)d + cj)\P) 

= —(P+-P+) 
M A { s ' 

P+ 

= 1 — ~r§~ — 1 (D.26) 

M A v ; 

where T ++ is one component of the canonical energy momentum tensor, P s + is the plus 
momentum of the scalar meson fields, and P + is the total nuclear plus momentum. The 
result Eq. (|D.26|) constitutes an analog of the Hugenholtz-van Hove theorem for finite 
relativistic systems; the equality becomes exact in the nuclear matter limit, where the scalar 
meson contribution vanishes, as shown in our previous work $9 . This means that we may 
anticipate that the binding effect will again be small. The vector operator 'mixing' and the 
scalar meson contribution will be elaborated on in a more general context in Section [D. 51 

It is also worthwhile to explicitly evaluate the expression Eq. (ITT231) for f(y) in the limit 
of infinite nuclear matter. In this case, V° = V + = V~ are constant and V± = 0, so we find 

/•oo 

A(z,x ± ) = / dz'V (z',x ± ) 

J z 

= -V°z 

= -V + z, (D.27) 

so that Eq. (iDTll becomes 

^+(fc) = J dxe' lk ^e- iik+ ~^ v+)z ^>+{x) 

= H(k±,k + -g v V+). (D.28) 

Therefore Eq. (lD~23ll b ecomes 

f(y) = 2M N ^2 f d 2 x ± \y+(x ±1 M N y-g v V + )\ 2 , (D.29) 
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which is simply the expression (|D.20|) modified by a shift in the argument of g v V + / Mm ■ 
Thus we find Eq. ()D.21|> is satisfied in the nuclear matter limit. It is important to stress that 
all that is recovered here is the shift in the argument and not any particular form of the plus 
momentum distribution which arises from the specific model used for the wavefunctions. 

D.3.2 Nuclear Matter 

In nuclear matter, we can derive an explicit form of f(y) valid for all mean field theories 
with fcp and Mn as the only input parameters. The manifestly covariant form of the single 
nucleon Green's function has been known for a long time |95| . and its use (in the nucleus 
rest frame) leads to the result 



X (k,P) = -i(j-(k-g v V)+M^) 

ITT 

(P-30) 



' + ^FT-A k ° ~ E *( k ) ~ 9vV°))6(k F - |k|) 



(k — V) 2 — Ml] + it E*(k) 



where 



E*{k) = JM* N 2 + k 2 . (D.31) 



The general form of the Green's function depends on a vector potential V = (V°,0) for 
a nucleus at rest, and the effective mass Mjy which includes the effects of interactions on 
the nucleon mass. The values of V and Mjy depend on the specific Lagrangian employed, 
but the form of the Green's function is general. Recall also that V~ = V + = V° for the 
expectation values of vector meson fields in the nucleus rest frame. 

The result Eq. (|D.30|) was first obtained using the conventional equal time approach, 
but the very same can also be obtained from the light front formalism. In that case it is 
necessary to include the effects of the instantaneous part of the nucleon Green's function 
and those of the instantaneous meson exchange. 

The next step is to insert the connected part (second term) of Eq. (|D.30|) into Eq. (|1.12|) 
for f(y). This gives, after taking the trace and using the delta function to integrate over 
k°, the result 

ff \ 4 f*u ^ E *( k ) + k \ , E*(k)+g v V + + k\ Q(h 

fiy) = Wfp~ B J dk ^ dk E*( k ) s{y w N )9{kF ~ |k|) - (D - 32) 
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The integration is simplified by using the transformation 

k + = E*{k) + k 3 , (D.33) 

which makes a connection with light front variables. It is an exercise in geometry to show 
that the Fermi volume can be re-expressed in terms of k + using 



k 2 ± + (k + -E*(k ± ,k + )) 2 <k%, E* F = y 'k% + M* N 2 , (D.34) 

k 2 + M* 2 ' 



(D.35) 



so that Eq. (|D.32|) becomes 

f{y) = (fcffe / ^ J dk+ ° ^ ~ kl ~ {k+ " E * {k))2) 5{V " k++ M^ V+ ) - (D ' 36) 
The use of the definition of the energy of a nucleon at the Fermi surface, 

E F = E* F + g v V + =E* F + g v V°, (D.37) 



allows one to achieve a simple expression for f(y): 

3 M 3 - r i.2 

f(y) = l^f e ^ E F + k F )/M N - y)8(y - (E F - k F )/M N )) 



(D.38) 



Mn m n 

The essence of the relativistic version of the Hugenholtz-van Hove theorem |62j follows 
simply from the statement that the pressure of a stable nucleus vanishes, and hence P + = 
P . Therefore, in the nuclear rest frame 

E A {k F ) = i(p+ + p-) =P+ = p- =M A 

E A (k F ) Ma 
A A 

E(k F ) = E F = M N (D.39) 
The result Eq. (|D.38|) can be further simplified to 

f(y) = -^e(l + A F - y)8(y -1 + A F ) [A 2 F - (1 - y) 2 ] , (D.40) 
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where Ap = Jzf/Mn, correct for any relativistic mean field theory of infinite nuclear mat- 
ter. Different theories with the same binding energy and Fermi momentum may have very 
different scalar and vector potentials, but must have the same f(y). 

A result very similar to Eq. (|D.40|) was previously obtained by Birse |13| . The difference 
between his formula and ours is the appearance of Mjv in the function f(y), whereas he 
uses Mat. This difference is a small effect numerically, and therefore our conclusions will be 
the same as his. 

The baryon sum rule and momentum sum rules are derived by taking the first two 
moments of f(y). This gives: 



The latter equation is remarkable; it states that in deep inelastic scattering the nucleons 
act as if they carry all of the P + of the nucleus even though the mesonic fields may be very 
prominent. 

This is clearer if we re- interpret these sum rules in terms of a probability f(k + ) that a 
nucleon has a plus momentum k + = yMjv> with f(k + ) = Af(yM]y)/MN, so that 



The momentum sum rule Eq. QD.44JI shows the total plus momentum carried by the nucleons 
(as seen in deep inelastic scattering) is also the total momentum carried by the nucleus. 
The main result of this is that the nuclear structure function is given by 11 . 12ft with 

the function f(y) obtained in ID.40|) . This tells us that, despite the fact that there is 

(A) 

considerable binding energy, there is no EMC binding effect. Indeed, -F 2 depends on the 
Fermi momentum but does not depend on the effective mass Mjy-. 




(D.41) 



(D.42) 




(D.43) 



(D.44) 
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D.4 Nuclear Structure Functions 

We determine the wave functions appearing in Eq. (|D.23|) numerically from a relativistic 
self-consistent treatment following Horowitz and Serot pl^ using the same program jOH] 
which includes electromagnetic effects. The plus momentum distribution follows and is 
given in Fig. ID. 31 for 16 0, 40 Ca, 208 Pb and in the nuclear matter limit. One can see that 
the peaks appear near y = 1 as required by the Hugenholtz-van Hove theorem Eq. (jD,26|) , 

It is worth noting that application of the Soper-Yan transformation Eq. (|D.3|) to the 
wavefunctions obtained from the equal time wavefunctions reproduces the plus momentum 
distributions, including the correct asymmetry, of the light front calculations in Ref. |21j . 
which did not use the approximation Eq. (jD.18|) . We show the comparison of our Oxygen 
calculation to that of Ref. |2J in Fig. ID.41 the agreement of these two curves demonstrates 
the excellence of the approximation relating light front and equal time wavefunctions. One 
can see that the effect in finite nuclei of the Soper-Yan transformation is to shift and broaden 
the plus momentum distribution, while in nuclear matter (also shown in Fig. ID.4)) it is just 
a shift. If these distributions were to be used in the nuclear structure function Eq. (|4.9|) 
though, since (y) ~ 0.8 for Oxygen, the ratio function (Eq. (|4.8|) discussed later) would fall 
precipitously to nearly zero at x ~ 0.6 in stark contradiction with experiment. 

The structure function is given by the convolution Eq. 1)4.9(1 with the generic function 
q = i*2 ^ • The assumption that nuclear effects do not modify the structure of the nucleon 
is embodied in Eq. ()4.9|) by the use of the structure function of a free nucleon; we use the 
parameterization |33| 

F 2 (A ° (x) = 0.58^(1 - x) 2 8 + 0.33 y/x(l - x) 3 ' 8 + 0.49(1 - x) 8 . (D.45) 

The experiments measure the ratio function Eq. 14.81 The results of our calculations are 
plotted for 16 0, 40 Ca, 208 Pb and in the nuclear matter limit in Fig. ID. 51 showing data for 
Carbon, Calcium and Gold from SLAC-E139 52 and an extrapolation |97| for the nuclear 
matter calculation. The most striking result is that these calculations fail to reproduce the 
EMC effect; the curves consistently miss the minima in the data, and the agreement gets 
worse with increasing A. Another important result is that the ratio function does not fall 
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Figure D.3: Plus momentum distributions, f(y), for 40 Ca (solid), 16 (short dashes), 208 Pb 
(dot-dashes) and nuclear matter (long dashes). 




Figure D.4: f(y) distribution for 16 (solid) after application of the Soper-Yan transfor- 
mation along with the 16 O distribution from Ref. j^J (short dashes). Note that the peaks 
occur at y < 1. 
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Figure D.5: Ratio functions for 16 0, 40 Ca and 208 Pb showing data for Carbon, Calcium and 
Gold, respectively, from SLAC-E139 52 . The nuclear matter calculation shows extrapo- 
lated data [97j . 
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to zero as would be the case if the small effective mass (~ 0.56Mjv for nuclear matter in 
QHD-I) were the relevant parameter describing the binding effect which would follow from 
using Eq. (|D.20|) instead of Eq. (|D.23|) . The results also show a minimum near x ~ 0.6 for 
Oxygen and nuclear matter that is deeper than the Calcium and Lead calculations. This is 
a curious feature that contradicts the trend in experimental data, and is due to the effects 
of two parameters. 

The first, and most important, is that of the location of the peak of the plus momentum 
distribution given by Eq. (|D.26|) . which gradually approaches y = 1 as the nuclear matter 
limit is reached. This is due to the fact that scalar mesons carry a small amount of plus 
momentum |21| that vanishes as A —* oo. The closer to y = 1 the peak is in Fig. lD.3| the less 
pronounced the minimum in Fig. ID.51 all else remaining constant. The second effect is due 
to Mjv, which reaches a minimum at 56 Fe corresponding to a more pronounced minimum 
of the ratio function than for A < 56 or A > 56, keeping the scalar meson contribution 
constant. 

Using a Gaussian parameterization of the plus momentum distribution and the experi- 
mental binding energy per nucleon via the semi-empirical mass formula, we have modeled 
the dependence of the minimum of the ratio function, R(x ~ 0.72), on the number of nu- 
cleons in the nucleus in Fig. ID.61 The motivation for the use of Gaussian plus momentum 
distributions is based on the expansion ^B] 

(N) 



Fi A \ XA ) = Ff\ XA ) + ex A dF - 



Ox 



X=XA 



+7 



2x A 



0f!; N) 2 2 F^ N) 

+ x A 2 

x=x A OX 1 



Ox 



X = Xj\ 



(D.46) 



where e = 1 — J dyyf(y) (D.47) 
1=1 dy(y - l) 2 f(y) (D.48) 



The Gaussian parameterization uses the peak location and width, (y) and ({y 2 ) — (y) 2 ) 1 ^ 2 
respectively, from the relativistic Hartree calculations in Fig. lD.3( and is normalized to unity. 
This allows us to obtain a plus momentum distribution for any A with minimal effort. We 
show the combined effect of scalar mesons and binding energy per nucleon on the ratio 
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Figure D.6: R(x = 0.72) as a function of A including scalar meson and binding effects (solid 
line), and leaving binding energy per nucleon constant at —8.5 MeV (dashed line). The 
data are from SLAC-E139 |H2]- 
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function along with the effect of scalar mesons alone using a constant binding energy per 
nucleon of —8.5 MeV independent of A. It can be seen that a changing Mn with A has the 
most effect for nuclei much larger than Iron, but does not change the general trend that the 
minimum of the ratio function becomes less pronounced as A increases due to the vanishing 
scalar meson contribution and the peak of the plus momentum distribution approaching 
unity. This dependence of the binding effect on A is quite different, both in magnitude 
and shape, than the trend in experimental data summarized in Ref. [HZj which satisfies 
R(x ~ 0.72) ~ A- 1 / 3 , so that the minimum becomes more pronounced as A increases. This 
fully demonstrates the inadequacy of conventional nucleon-meson dynamics to explain the 
EMC effect. 

D.5 Scalar Mesons and More General Considerations 

The average value of y, given by Eq. (|D.26|) . yields the nucleon contribution to the plus 
momentum, and is less than one which can be seen in Fig. ID.3I We now address the 
remaining plus momentum in finite nuclei. Previous results [2^ show that a small fraction 
(5y ~ 0.005) of the plus momentum is carried by the scalar mesons which vanishes as the 
nuclear matter limit is approached. This is due to the fact that scalar mesons couple to 
gradients in the scalar density (arising mainly from the surface of finite nuclei) which vanish 
as A — > oo. The question is: why are scalar mesons allowed to carry plus momentum and 
not vector mesons? 

The simplest answer lies in the Dirac structure of Eq. 1)1.12(1 : the 7 + in the trace picks 
out terms in the full interacting Green's function with an odd number of gamma matrices 
which includes all Lorentz vector interactions and excludes Lorentz scalar interactions. The 
Dirac structure of f(y) is directly related to the Dirac structure of the energy momentum 
tensor, so the answer also lies there and illuminates a problem with conventional nucleon- 
meson dynamics. The component of the energy momentum tensor relevant to the plus 
momentum, from a chiral Lagrangian containing isoscalar vector mesons, scalar mesons and 
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pions, is given by [ZSl EH] 

+ 9+7T • + 7T • d+TT _ - (1 - ^ sill 2 ~). (D.49) 

Since each of the terms in Eq. ()D.49|) involves one of the fields, it is natural to associate each 
term with a particular contribution to the plus momentum. This decomposition, though, is 
not well defined; field equations relate various components. We see the first three terms of 
Eq. I)D.49|) appear in (y), which defines the nucleon contribution to the total nuclear plus 
momentum, in the derivation of the Hugenholtz-van Hove theorem Eq. (|D.26|) : we are not 
allowed to have the vector mesons contribute a well defined fraction of plus momentum. This 
means that one could trade certain mesonic degrees of freedom for nucleons by replacing 
mesonic vertices with nucleon point couplings, for example, in line with the general concept 
of effective field theory. In our case the first three terms are related by the vector meson 
field equation, but the fourth is left out since the scalar mesons couple to the scalar density 
*I n I' which is not present in Eq. (|D.49|) . Therefore the scalar mesons (and pions) contribute 
a well defined fraction of plus momentum. These explicit meson contributions create an 
EMC binding effect, but the pionic contributions are also limited by nuclear Drell-Yan 
experiments [2] to carrying about 2% of the plus momentum which is insufficient to account 
for the entire EMC effect which corresponds to about 5% of the plus momentum for Iron. 
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